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An operator convex function on (0, oo) which satisfies the symmetry condition k(x^ 1 ) = 
xk(x) can be used to define a type of non-commutative multiplication by a positive definite 
matrix (or its inverse) using the primitive concepts of left and right multiplication and 
the functional calculus. The operators for the inverse can be used to define quadratic 
forms associated with Riemannian metrics which contract under the action of completely 
positive trace-preserving maps. 
. — We study the question of when these operators define maps which are also completely 

^ positive (CP). Although A ^ D^ 1 / 2 AD^ 1 / 2 is the only case for which both the map and 

its inverse are CP, there are several well-known one parameter families for which either 
the map or its inverse is CP. We present a complete analysis of the behavior of these 
families, as well as the behavior of lines connecting an extreme point with the smallest 
one and some results for geometric bridges between these points. 

Our primary tool is an order relation based on the concept of positive definite func- 
tions. Although some results can be obtained from known properties, we also prove 
new results based on the positivity of the Fourier transforms of certain functions. Con- 
crete computations of certain Fourier transforms not only yield new examples of positive 
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definite functions, but also examples in the much stronger class of infinitely divisible 
functions. 
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1 Introduction 

On a commutative algebra the operations of multiplication and "division" by elements of 
the positive cone take the positive cone into itself. However, this is not the case for non- 
commutative algebras, on which these operations are not even uniquely defined. 

Various non-commutative versions of multiplication and division (i.e., multiplication by 
the inverse) by elements of the cone of positive definite matrices or operators correspond to 
maps on matrix algebras, and some of these maps play an important role in many contexts. 
For D > 0, some naive definitions of multiplication by the inverse are given by the maps 
X H> D^^XD- 1 / 2 , X ^ D~ l X and X H- ID" 1 . The first maps the cone of positive 
operators to itself, while the other two do not even preserve self-adjointness. There are many 
other possible generalizations of multiplication by D^ 1 . In Section [Z] we consider several 
different one parameter families of such maps. 

Perhaps, the best known example is 



D 



(A) = / (D + tI)- l X(D + tI)~ l dt. (1.1) 
Jo 



which gives a well-defined (and highly symmetric) notion of non-commutative multiplication 
by D^ 1 . Its inverse is well-known to be 

n~ 1 (Y)= [ D t YD 1 ~ t dt (1.2) 
Jo 

which is a form of non-commutative multiplication by D. The quadratic form Tr A*£l~j^~(B) 

is known as the Bogoliubov or Kubo-Mori inner product. Although the inverse relationship 



(1.1) between (1.2) is well-known, and follows from more general results given later, most 



proofs in the literature, e.g., |344 Section 3] rely on expansions in eigenvectors which are not 



very satisfying. Therefore, we include an explicit proof in Appendix C.2 

We consider here a class of such maps which arise in quantum information theory, in the 
context of what are known as monotone Riemannian metrics |45|, [3"3"] . 

We study the question of which maps within this class have the property known as com- 



pletely positivity defined in Section 2.3 The map in (1.1) has this property, but its inverse 



(1.2) does not. 
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This question was motivated by an observation in |33j E] about bounds on the contrac- 
tion of monotone metrics under the action of completely positive trace preserving (CPT) 
maps which are also known in quantum channels in view of their important role as noise 



models in quantum information theory. These bounds are discussed briefly in Section 2.4 
and Appendix [B| 

This paper is organized as follows. Section [2] describes the various concepts we need 
and introduces the notations we will use. Subsection |2.4| contains a brief summary of the 
background and motivation behind this work. In Section [3] we present some powerful tools we 
will use formulated in terms of positive kernels, and a closely related partial order. In Section|4] 
we present a large number of examples of one parameter families of functions which provide 
a number of inequivalent classes of maps used to define non-commutative multiplication by 
the inverse of a positive matrix. In most cases, we can also provide precise ranges for which 
these maps or their inverses are CP. Section [5] proves new results about positive kernels based 
on Fourier transforms, which are needed for some of our results. These results are of interest 
in their own right. In Section [6] we complete the proofs of those results stated in Section [4] 
which require the results of Section [5] 

There are also three appendices. The first contains the detailed proof of the integral 
representation needed in Section 2.1 The second gives more details about the motivation in 
terms of contraction under CPT maps of the Riemannian metrics associated with our maps. 
Finally, for the benefit of non-experts, we present some very pedestrian arguments which 
clarify some well-known results that are often glossed over. 



2 Preliminaries 

2.1 Basics 

For each d E N we write M<2, H<2, and P^ for the sets of d x d complex, Hermitian, positive 
definite and positive semi-definite matrices, respectively. Functions of matrices in can be 
defined by using the spectral theorem; this is sometimes called "functional calculus" (see, 
e.g., 091 Section VII. 1]). 

A real function / on (0, oo) is said to be operator monotone (or operator monotone in- 
creasing) if A > B implies f(A) > f(B) for every A,B E P^ with any d E N, operator 
monotone decreasing if — / is operator monotone. A real function k on (0, oo) is said to be 
operator convex if 

k(XA + (1 - X)B) < Xk{A) + (1 - X)k(B) 

for all A, B E Pd with any d E N and all A E (0, 1), and operator concave if —k is operator 
convex. The theory of operator monotone and operator convex functions was initiated by 
Lowner and Kraus [23] , respectively. It is well-known Section V.4] (also [HHHE!]) that 
operator monotone (also operator convex) functions on (0, oo) have an analytic continuation 
into the upper half-plane of C. Moreover, a necessary and sufficient condition for a function 
on (0, oo) to be operator monotone increasing (resp., decreasing) is that it has the "Pick" 
mapping property that the analytic continuation maps the upper half-plane into the upper 
(resp., lower) half-plane. The integral representation theory for Pick functions and operator 
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monotone functions are also well-known, see, e.g., [U Section 59] or [TJ Section V.4] and 

[21E21HH]. 

Definition 2.1. Let K, denote the class of functions k : (0,oo) — > (0,oo) which are operator 
convex and satisfy the symmetry condition xk{x) = /c(x _1 ) and the normalization condition 
fc(l) = l. 

There are a number of equivalent characterizations of the class tC which are given in 



Theorem 2.4 below. Its proof uses an integral representation, which is important in its own 



right and presented in Theorem 2.3 We also observe that although k(x) 6 K, may diverge as 
x \t 0, it cannot diverge more rapidly than x -1 . This was proved in [22] • For completeness, 
we include its proof as well as the proof of Theorem 2.3 in Appendix [A) 

Proposition 2.2. Let k : (0, oo) —> (0, oo) be operator convex. Then lim^^o xfc(x) exists and 
is finite. When xk(x) = £;(x _1 ), lim^-^oo k(x) also exists and is finite. 

Theorem 2.3. For every k & K, there exists a unique probability measure m on [0, 1] such 
that for x £ (0, oo) 

1+x (1 + ^) 2 



k(x) = [ r • - — — — — dm(v 

J[o,l] (x + v)(l + vx) 2 



[0,1] 



1 + r^^^M. 



x + v 1 + xv ) 2 
Theorem 2.4. For each function k : (0, oo) —> (0, oo) the following are equivalent: 

(a) k is operator convex and xk(x) = fc(x _1 ); 

(b) k is operator monotone decreasing and xk{x) = fc(x _1 ); 

(c) f{x) = l/k(x) is operator concave and f(x) = x/(x _1 ); 



xfix r ) 



(d) f(x) = l/k(x) is operator monotone and f(x) = xf( 

Proof. The equivalence of xk(x) = fc(x _1 ) and f(x) = x/(x~ 1 ) is easily checked. The 
equivalence (c) 43- (d) for positive functions on (0,oo) is well-known (see, e.g., [TJ Theorem 
V.2.5]), and (b) 44> (d) follows from the fact that x i— > 1/x is operator monotone decreasing on 



(0, oo). The implication (a) (b) follows immediately from Theorem 2.3 and the well-known 
fact that the map x h-> l/(ax + f3) is operator monotone decreasing on (0,oo) for any fixed 
a,/3 > 0. We finally show that (b) =^ (a). Assume (b); then l/k(x) is operator monotone and 
so operator concave on (0, oo). This implies (a) since x _1 is operator monotone decreasing 
and operator convex. QED 



As shown in the above proof, the implication (b) (a) holds true without the symmetry 
assumption xk(x) = ^(x^ 1 ). However, the reverse implication (a) =^ (b) only holds under 
this additional assumption and appears to be new. 

The next result is easy to verify, but stated explicitly for completeness. 

Proposition 2.5. The map k i-> k = 1/A;(x _1 ) is a bijection on K, and the map f i— > f = 
l//(x _1 ) induces the same bijection with f(x) = l/fc(x) and /(x) = l/k(x). 
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2.2 The multiplication map and its inverse VL k D 

For each D 6 we write Ld and Rd for the left and the right multiplication operators, 
respectively, i.e., LpX = DX and RdX = XD for X G M^. Note that Ld and Rd are 
commuting positive invertible operators on considered as a Hilbert space equipped with 
the Hilbert-Schmidt inner product (X, y)ns — TiX*Y, where Tr denotes the usual trace 
functional on M^. The operator LaR^ 1 was used by Araki [5J to define the relative entropy 
of positive operators A, B in far more general situations than matrix algebras, and is often 
called the relative modular operator. 

For a fixed function k G K, we define, for any D G P^, the linear map 0^ : — > by 

n k D (X) = R-Jk (L D R D l ) X = L D x k {RdLo 1 ) X, X E M d . (2.1) 

Since both of the commuting operators Ld and Rd are positive with respect to the Hilbert- 
Schmidt inner product, it is clear that fl D is also positive (in the same sense). Each map 
Sl D can be considered as a non-commutative generalization of multiplication by D -1 ; indeed, 
if DX = XD then Q D (X) = D~ 1 X independently of k G K. The equivalence of the two 



expressions in (2.1) follows from the symmetry condition xk(x) = k(x x ) since 
R-Jk {LdR-J) = L^LoR^k (LdRb 1 ) 

= L n l k ((LdR-^V 1 ) = L^k (RdLJ 1 ' 



To better understand the action of Q D , we consider the two- variable function 
(p k (x,y) = {1 / y) k{x / y) for x,y > and observe that Q, D = 4> k (L d , Rd) ■ When D is a 
diagonal matrix with eigenvalues \j, it is an easy consequence of the functional calculus that 
the action of Q D on the matrix with entries is 



4> (Xj, \k)xij = —k 



1 , / A 



a j V Aj 

which is the Schur (or Hadamard or pointwise) product AoX with the matrix A with entries 
a>ij = fc (Aj, Xj). More generally, let U be a unitary which diagonalizes D so that 

D = Udiag(X 1 ,...,X d )U*. 

Then 

n k D (x) = u([4> k (\i, \j)} o [u*xu])u*. (2.2) 



Since L^, 1 = L^-i and R D ^ = Rd- 1 , it might be tempting to think that ^^(X) = 
rj£)-i(X). However, this is easily seen to be false by considering specific examples (including 



(1.1) and (1.2)). Instead we have for any D G IPrf, 

J f D = (^d)' 1 = RDf{L D R D l ) = RdHL^Rd) = n % D ^ (2.3) 

with f(x) = l/k(x) = /c(x _1 ). To see this, observe that it follows from (2.2) that (il^,) -1 = 
(l/4> k )(LD,RD)- From the relation 

1 ..s _ V 

k(x/y) 



uk ( x ' v) = uf-i..\ = yf( x /y) = x f(y/ x )> 



the functional calculus implies (2.3) 
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2.3 Complete positivity of 

A linear map <3? : i— > is called positive if it is positivity-preserving in the sense that 
A > implies > 0, i.e., <b(F d ) C P d . A linear map $ : M d i-> M d is called completely 

positive (CP) if ® X„ is positive on £g> M n for all n G N with X n the identity map on 
M n . The notion of complete positivity, introduced by Stinespring [53J and discussed in, e.g., 
|42l Chapter 6] plays an important role in quantum information theory. (See, e.g, \AQ \ 147 ) .) 

The recognition in (2.2) that ttj-, can be implemented as a Schur product yields a simply 
stated condition to test that it is CP. In general, complete positivity of a map on is a much 
stronger condition than positivity. However, for Schur products, it is well-known (see, e.g., 
[4"2l Theorem 3.7]) that both positivity conditions for the map = A o X (A, X G M d ) 

are equivalent to positivity of A. Indeed, the map $a ® 

X n on M n (Md) = Md ® M n can be 
realized as Schur multiplication with A<g> J n , where J n is the n x n matrix with all entries 1. 
Therefore, in our setting, the requirement that the map is CP is equivalent to the weaker 
positivity requirement, as we explicitly state for completeness in the following: 

Theorem 2.6. The following conditions for k G K, are equivalent: 

(a) : M d -> M d is CP for every D £ F d with any d G N; 

(b) fi^J : Md — > Md is positive for every D G Fd with any d G N; 

(c) the d x d matrix 

A ■ 



(2.4) 



l<i,j<d 



is positive semi-definite for every wi, . . . , > with any d G N. 



The next result allows us to replace the matrix A in part (c) of Theorem 2.6 by some 
closely related matrices for which the positivity condition may be more easily checked in 
some situations. 



Proposition 2.7. The matrix A in (2.4) is positive if and only if one (and hence both) of 
the following matrices are positive: 



Wjk 



Wi 



l<i,j<d 




Wi I Wi 

k 



(2.5) 



l<i,j<d 



Proof. Let W be the diagonal matrix with entries WiSij. Then the matrices above correspond 
to W*AW and (W*) 1 ' 2 AW 1 / 2 , respectively. QED 

Indeed, note that the condition xk(x) = ^(x^ 1 ) implies that A = W* AW . 



2.4 Background and motivation 

For each k G JC, the map Q k D can be used to define a quadratic form 

T k D (X,Y) = (X,Q k D (Y))ns = TrX*Q k D (Y) 



(2.6) 
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which can be interpreted as a metric on the Riemannian manifold T> d = {f G : Tr D = 1} 
of invertible density matrices in M. d . Here, the matrices in M d with trace zero form the 
tangent space, denoted by M d , of V d at each foot point D. This metric is monotone in the 
sense that for any completely positive and trace-preserving map $ : — > M m (d, m G N), 

r| (D) ($pQ, cD(y)) < r k D (x, y), d g v d , x, y g m°. (2.7) 

The theory of monotone Riemannian metrics was largely developed by Petz |45j after Moro- 
zova and Chentsov [3H] introduced the concept. It was shown in [33] that each k G fC defines 



a family of monotone metrics of the form (2.6) with D 6 T> d for all d G N, and that any 



Riemannian metric on T> d , d € N, which satisfies the contraction condition (2.7) must be of 



the form (2.6) for some k G /C. (See also [32].) 

In |45| . the operator Jjrj defined in (2.3) was used to define monotone metrics in the 



equivalent form as 

r k D (x,Y) = (x, (4)- 1 (y)) HS = TTX*(jf D y\Y). 

It might seem more natural to work with J D which is a non-commutative version of multi- 
plication by D rather than using its inverse (introduced in [33]). However, in this paper 
we use Q k D instead of jj^ since it avoids the need to take inverses to define our target maps. 

In |54j . monotone metrics of the form Tq(P— Q, P—Q) with P, Q G T> d played an important 
role in the study of mixing times of Markov processes. It was observed in j54l Section III.B] 
and [33l Section IV.C] that when both £l k D and its inverse (O^) -1 are positivity-preserving, 
one can obtain a useful upper bound on the contraction of Riemannian metrics, which is 
described in more detail in Appendix pi In |54] this bound was used in the case k(x) = x -1 / 2 
for which both n k D (A) = D^^AD^/ 2 and the inverse D l / 2 AD 1 / 2 clearly map Prf into 



itself. (In fact, for every D G IP^ they are bijections on P^-) Theorem 3.5 below implies that 
k(x) = x -1 / 2 is the only function in K, with this property. 

The study of quasi- entropies was also initiated by Petz in [33l HH HI] , which can be defined 
from any operator convex function g on (0, oo) with g(l) = as 

H g (A,B,K) = (K^^aR-^RbK)^ = Tr \* r BK* g^LAR^ 1 ) (KVB) (2.8) 

for A, B G ¥ d , K G M d . It was later observed in [HH 0H] that for any D G V d the Hessian 

d 2 



dadb 



H g (D + aX,D + bY,I) 



a=b=0 



X,Y£E° d , 



can be associated with a monotone Riemannian metric is some important examples. This 
was proved for more general g in [33, Theorem II. 8], where it was also noticed that g(x) = 
(1 — x) 2 k(x) with k G K, is an operator convex function on (0, oo) with <?(1) = 0. Moreover, 
the symmetry condition xk(x) = /c(x -1 ) implies that xg(x _1 ) = g(x) and the quasi-entropy 
with K = I has the symmetry property 

H g (A,B,I)=H g (B,A,I). (2.9) 

and that every quasi-entropy with this symmetry property comes from a k G fC. The quantity 
H g (A, B,I) is often called an f -divergence. See [21] for a thorough discussion of /-divergences 
(without the additional symmetry condition). 
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When k(x) = 4/(1 + ^fx) 2 so that g{x) = 4(1 - ^/x) 2 , the function H g (A,A,K) is the 
Wigner- Yanase skew information [56] , which Dyson suggested extending to the case including 
the parameter p £ (0, 1), that is equivalent to using g(x) = 4(1 — x p )(l — x 1 ~ p ). This led 
to Lieb's seminal work on concave trace functions [33], in which he showed that (A,B) \— )■ 
Tr K*A P KB 1 ~ P is jointly concave in A, B £ whenp G (0, 1). It is implichj^jin Ando's paper 
[3J that the quasi-entropy H g (A, B, K) can be extended to g(x) = (1 — x p )(l — x l ~ p )/p(l —p) 
with p £ [—1, 2]. Hasegawa |17) seems to have been the first to use well-known properties of 
monotone and convex operator functions to explicitly recognize that replacing 4 by l/p(l —p) 
allows one to extend the quasi-entrop}|^] for the WYD skew information and the associated 
Riemannian metric to the maximal range p 6 [— 1, 2] (with p = 0, 1 defined as limit^. See 
also [23] where equality conditions were given for the convexity of H g (A, B, K) and some 
other inequalities for the extended WYD family. 

In this paper, we make use of tools developed by Hiai and Kosaki [19, 20J in study of 
means of operators]^] Motivated by this work, whenever k S /C, we define 

Mk{x > y) ^WJ7y x ' y>0 ' (2 - 10) 

M k (A,B) = R B (k(L A R B 1 )Y\ A,BeP d . (2.11) 



From (2.1) and (2.3) we have in particular 

M k (D, D) = (Sl'hy 1 = fi* = J f D . (2.12) 

The function M(x,y) is called a symmetric homogeneous mean for positive scalars, i.e., M : 
(0, oo) x (0, oo) — > (0, oo) is a continuous function such that 

(1) M(x,y)=M(y,x), 

(2) M(tx,ty) = tM(x,y) for t > 0, 

(3) M(x, y) is non-decreasing in x, y, 

(4) min{x,y} < M(x,y) < m&x{x,y}. 



With f = 1/k, definition (2.10) is equivalent to M (x, y) = y f(x/y) which was used in 



under the weaker condition that / is non-decreasing in the numerical sense. It follows from 



Proposition 2.5 that as k runs through K, both conventions generate the same set of operators 



of the form (2.11) 



1 Ando found an alternate proof of Lieb's concavity results and also showed convexity for p S (1, 2). Both 
Lieb and Ando ignored the linear term TrK*AK in the skew information, since it was irrelevant to convexity. 
2 Hasegawa actually used the asymmetric g(x) = (1 — x p )/p(l — p). However, it follows from Eq. (37) in 



[33| that this yields the same k(x) given by (4.61 as the symmetric version above. 

3 Lindblad [35i was the first to observe that one could recover joint convexity of the usual relative entropy 
by taking limp_>i in Lieb's result. 

4 This work was motivated by inequalities for unitarily invariant norms. The term mean used there does 
not, in general, yield the mean of a pair of operators in the sense of Kubo and Ando [31] . 
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2.5 The convex sets K, and /C + 

Recall that K, denotes t he s et of functions k : (0, oo) —> (0, oo) satisfying any of the equiv alent 



2.4 



and k(l) = 1. With k(x) = l/k( X - 1 ) given in Proposition 



2.5 



conditions of Theorem 
k \— > k is a bijective transformation on K. 

For i/ G [0, 1] let us set 

y^Mig. \+- -fl±iO('-J- + ^_y (2.13) 
2 (a; + i/)(l + i/a:) 2 Vx + i/ 1 + 2*7 V ; 

For any fixed x G (0, oo), by computing the derivative of k v {x) in v one can easily verify that 
k v {x) is non-increasing in v G [0, 1] so that 

= y^< K x \x) < = *T * e (o, l). (2.14) 

Since Theorem 12.31 can be rewritten as 

k( x ) = [ k cxt {x) dm{v) (2.15) 
J[0,1] 

with m a probability measure, one moreover has 

< k(x) < k£lC. (2.16) 



1 + x 2x 

Thus, K, has the smallest element fc| xt (x) = 2/(1 + x) and the largest element A;Q Xt (x) = 
(1 + x)/2x in the pointwise order. 

Now we may consider K, as a subset of the locally convex topological vector space consisting 
of real functions on (0, oo) with the pointwise convergence topology. Then it is obvious from 



(2.16) that /C is a convex and compact subset. The uniqueness of the representing measure m 



in Theorem 2.3 implies that /C is a Choquet simplex with the extreme points kf, Kt for v G [0, 1] 
(that is the reason for the notation k° xt ). Furthermore, since v i— > kf, xt is a homeomorphism 
from the interval [0, 1] into /C, one sees that K, is a so-called Bauer simplex (as in |16j). 

Motivated by the work on contraction bounds in [331 El] which is described in Appendix [fH 
we define two subsets K, + and K~ of K, as 

JC + = {k G K : fl k D is CP for every DeF d ,d£ N}, 
KT = {k G K : (^o) _1 is CP for every D G P^, d G N}. 



In view of (2.9), an operator convex function k on (0, oo) satisfying xk(x) = k(x x ) might 



be called a "symmetric" operator convex function. It would then be natural to call functions 



in /C + "completely symmetric". It follows from (2.3) that 



k G /C + <J=^> k G /C~ where fc(x) = l/A;^ 1 ) . (2.17) 



It follows from from Theorem 2.6 that /C + and /C are closed under pointwise convergence. 



Although 1C + is convex, K, is not convex (as shown in Example 4.4 below). Since JC + is a 
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compact convex subset of IC, it is the closed convex hull of its extreme points by the Krein- 
Milman theorem. However, determining all the extreme points of IC + seems quite challenging. 



Some non-trivial ones are described in Example 4.3 

In this paper, we have chosen to formulate most of our results in terms of functions k £ IC. 
As is clear from Theorem |2.4| we can also define the convex set of functions F with f = 1/k 
which satisfy property (c) or (d). Although our choice is partly a matter of taste, in some 
situations, one may be more convenient than the other. We find it useful here to let 

F ± = {f £F: {n]i f ) ±l is CP for every D € Pd, d G N}, 

so that IC^ corresponds to F^ by k ■<->•/ = 1/k. Since l/k(x) = k(x _1 ) = x k(x), it is obvious 



by (2.17) that 



F = {xk(x) : k £ /C}, F + = {xk(x) : k £ KT}, F~ = {xk(x) : k £ /C+}. 

Hence k xk{x) gives an affine correspondence between fC and F, by which /C + is isomorphic 
to F~ . Therefore, F~ is also convex and the extreme points of F are 

3 Positive kernels and induced order 
3.1 Basic definitions 



In principle, the condition of Theorem 2.6 (c) gives a simple criterion for complete positivity. 



But in practice, it is not easy to verify that either the matrix A in (2.4) or one of those 



in (2.5) is positive semi-definite. Only a few examples can be resolved using this criterion. 
However, there is another equivalent condition based on the theory of functions which define 
positive kernels. 

Definition 3.1. A continuous function h : M i— > C is called positive definite if h(x — y) 
is a positive semi-definite kernel, i.e., [h(U — tj)] 1<i <d is positive semi-definite for any 
tl, . . . , td € K with any d € M, or equivalently, 



<p(s)h(s-t)(p(t)dsdt>0, e C °°(l), 

where Cq°(]R) denotes the smooth compactly supported functions on R. Functions satisfying 
this condition are sometimes called "functions of positive type" [501 Section IX. 2] or "positive 
in the sense of Bochner" . 

Moreover, h is called infinitely divisible if h{t) r is positive definite for every r > 0, or 
equivalently, /i(t) 1//n is positive definite for every n G N. 

For convenience, some basic properties of positive definite functions stated here: 

(a) A positive definite function h is uniformly bounded on M. as \f(t)\ < /(0). 
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(b) Bochner's theorem (see |50i Theorem IX. 9], [TJ Section 60]) says that h is positive 
definite if and only if it is the Fourier transform of a finite positive measure on M. 
Thus, positive definiteness of h can be checked, in principle, by testing positivity of its 
Fourier transform. 

(c) The product of positive definite functions is positive definite. This immediately follows 
from the well-known fact that the Fourier transform of the convolution of two finite 
measures is the product of their Fourier transforms, or from the Schur product theorem 
for positive semi-definite matrices. 

In this paper we only consider positive definite functions on M so that we shall omit "on 
W in the rest. Positive definite functions played an important role in the work [2"U] on 
means of operators, where a partial order was introduced. The following definition is its 
adaptation to functions in K,: 

Definition 3.2. For ki,k 2 G /C we write k\ =^ k 2 if either of the following equivalent condi- 
tions holds : 

(a) the function kiie 1 ) / k^e 1 ) is positive definite on M; 

(b) the matrix 



ki(wj/wj) 
k 2 (wi/wj) 



l<i,j<d 



is positive semi-definite for every w\, . . . , Wd > with any d G N. 

It is easily verified as in [19, 20J that ^ is really a partial order in /C, and k\ =^ k 2 implies 
k\{x) < k 2 (x) pointwise. 

The stronger condition that ki(e t )/k 2 (e ) is infinitely divisible (following Definition 3.1), 



was studied in [8j. Results given there sometimes play a role in showing that the one- 



parameter families studied in Section 4.2 are monotonic in the ==! order. Moreover, infinite 



divisibility is important in the discussion of geometric bridges in Sections 4.3 and 6.3 Some 
examples considered here require new results for specific functions which are obtained in 
Sections 15.31 and 15.41 



The next useful lemma on positive definite functions will often be used in this paper. See 
|25| Appendix B] and f9j Theorem 3.2] for the proof of (1). On the other hand, (2) was first 
proved in [91 Theorem 5.1] while the "if part" was pointed out earlier in |57| . 

Lemma 3.3. 

(1) The function sinhat/ sinht is positive definite for a G (0, 1). 

(2) For (3 > — 1, the function (coshi + /3) _1 is positive definite if and only if f3 < 1. 



3.2 Basic applications 



The next theorem gives a basic characterization of the class tC + . The equivalence of (a)-(c) 
follows immediately from Theorem 2.6 with A replaced by the second matrix in (2.5). The 
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equivalence of (b) and (d) is an adaptation of |19l Theorem 1.1] via (2.12) in the present 
situation. 



Theorem 3.4. The following conditions for k £ JC are equivalent: 

(a) k G JC + , i.e., is CP for every D G P^ with any d G N; 

(b) k^x- 1 / 2 ; 

(c) e t ^ 2 k(e t ) is positive definite; 

(d) i/iere exists a symmetric probability measure v on R snc/i i/iat 



(3.1) 



/or a// DeP n and X £M d with any d G N. 



It is a well-known consequence of the Stinespring representation theorem that a CP map 
$ on the matrix algebra can be represented in the form $(A) = Y2j Fj^F* with at 
most d 2 matrices Fj G M d (see, e.g., [301 [TO], [HI Proposition 4.7] or [Ml Appendix A]). 



Thus, for any fixed D G P<i, when 0-^ is CP, one can find matrices Fj in such that 
VL k D {X) = Yl'JLiFjXF? with m < d 2 . But, for fixed k G /C, the representation will change 
with L> (hence with d) . (Even for fixed D the in the representation are only determined up 
to a unitary transformation Fj i— )■ n ij-^ with ity entries of a unitary matrix.) However, 



we are allowed to use integral representation, Theorem 3.4 says that we have the standard 
representation given in (3.1), from which the CP of the map fl^ is directly seen. Moreover, 



one sometimes has different integral expressions of or (O^,) 1 ; a typical example is (1.1) 
for n k D in case of k(x) = logrr/(x — 1) (see Appendix C.2). 

It follows immediately from (2.17) and Theorem 3.4 that k G KL~ if and only if k )>= x" 1 ! 2 



Consequently, x 1 I 2 is the largest element of /C + and the smallest of K, 
is a partial order on JC, we conclude 



Moreover, since ^1 



Theorem 3.5. The only function in K, for which both VLp and O^, 1 are CP for every D G P^, 



de 



is x 



-1/2 



It follows from Theorem 3.4 that the the problem of determining whether or not k G JC 
belongs to JC + can be reduced to the computation of the Fourier transform of the function 
e t / 2 k{e t ). However, this is often a hard task as will be seen in Section [5] 

In contrast to /C, it does not seem easy to find extreme points of IC + other than x -1 / 2 
and 2/(1 + x) which are the largest and the smallest elements of IC + , respectively, in the 
order =<! as well as the pointwise order. However, some new extreme points will be described 
in Example 4.3 and Theorem |6.2[ In addition, a natural boundary point will be found in 



Example 4.8 which is conjectured to be an extreme point. 



By comparing part (b) of the next result to (2.14), one immediately sees that =^ is stronger 
than the pointwise order. 
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Proposition 3.6. The following relations hold. 
(a) kf xt (x) = T - rz k cxt and k cxt 4 ^J-^ = k^ xt (x) for all v G [0, 1]. 



(b) k ext 4 



1 + X 
1 + X 



2x 



2x 



kQ Xt (x) and kf xt (x) 



1 + x 



4 A£ xt for all v G (0,1). 



(c) 2/(1 + x) 4 x- 1 / 2 4(1 + x)/2x. 
Proof. A straightforward computation gives 

fcfV) _ A£ xt (e-*) 4 (e* + v)(l + z/e*) 

*£ xt (e*) ~ /fcg xt (e-') ~ (l + i/) 2 (1 + e*) 2 
4 j/(e* + e"< + 2) + (1 - v) 2 



(1 + 1/)! 



4i/ 2(l-i/) 



e* + e - * + 2 
2 1 



{l + u) 2 (l + i/) 2 coshf + 1 



from which (a) follows by using /3 = 1 in Lemma 3.3 (2). Similarly 

*£ X V) = fcf^Ce"') =( )2 e* 
/- xt (e*) jfe«t( e -*) 1 J (e* + i/)(l + i/e*) 
2 1 
- { - 1 + u > „( e t + e -t) + l + i/2 

_ (1 + ^) 2 1 

2i/ coshi + ±±^' 



Since (1 + z/ 2 )/2i/ > 1 for 1/ G (0,1), this proves (b) by Lemma |3.3| (2) again. Finally (c) 
follows easily from 



Jl 2 kf-\e l ) 



-t/2 



2e 1 / 2 



1 



fcgxt( e f) e t + i cosh(t/2)' 



QED 



Proposition 2.2 implies that for every k G JC, xk(x) is bounded on (0,6) and k(x) is 



bounded on (a, 00) for any a, b > 0. Theorem |3.4| implies that a necessary condition for 
G ZC + is the stronger property that x 1//2 /c(:r) is bounded on (0, 00). However, this is not a 
sufficient condition. Indeed, it holds for all kf, xt (x) with v G (0,1]. Yet, as will be seen in 
Example 4.1 k'^ t (x) G fC + only for v = 1. The following result will be used in Example 4.4 



to analyze convex combinations of x 1 / 2 and /c^ xt . 



Proposition 3.7. Assume that k G JC\ JC + and lmx^oo x 1 l 2 k(x) = 0. Then for every 
AG (0,1], 

Afc + (1- X)x- 1 / 2 <£ K+. 
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Proof. Assume that k G JC satisfies linx^oo x 1 ^ 2 k(x) = and Xk + (1 — A)x -1 / 2 G JC + with 
some A G (0, 1]. Then, thanks to Theorem 3.4 and Bochner's theorem there is a probability 
measure /i on i satisfying 

roo 

Ae*/ 2 A:(e t ) + (1 - A) = fi{t) = / e its d(x(s), t G K. 



However, the symmetry condition xk(x) = k{x l ) implies e t ^ 2 k(e t ) = e t / 2 k{e *), t £ 1, 
and hence hmw^oo e t ^ 2 k(e t ) = 0. Therefore, we have 

//({()}) = , |im £(i) = 1 — A 

|t|— >oo 

(see |20t Corollary A.8]). This means e t / 2 k(e t ) = p, Q (t), t G R, with the probability measure 
/io = A _1 (/i — fi({0})So), implying the contradiction k G JC + . QED 



4 Examples 

In this section we list known families of functions in JC and investigate which functions in 
those families belong to JC + (or JC~). In this way we will see that JC + indeed contains a 
variety of functions even though it occupies only a small part of JC. 



4.1 Extreme points and simple averages 

Example 4.1. (Extreme points of JC) The extreme points of JC are A£ xt , v G [0, 1], given in 
(2.13). These are not in JC + unless v = 1 for which we have A;f xt (x) = 2/(1 + x). Indeed, for 
v G (0, 1] we find 



(1 + v) 2 cosh(t/2) 
2i/ ' cosht+ 



If e t//2 /c£ xt (e*) is positive definite, then so is its product with the positive definite 1/ cosh(£/2). 
But this yields (up to a constant) a function of the form in Lemma 3.3 (2), which is not positive 
definite for f3 = (1 + v) 2 /2v > 1 when v G [0, 1). 

It was shown in [6., Example 9] that k^ it (x) < x -1 / 2 (in the pointwise order) for all x > 
if and only if 3 — 2\/2 < v < 1. This example provides another demonstration that the =^ 
order is stronger and =^ x -1 / 2 is the key to determining whether or not a function k G JC + . 

Example 4.2. (Convex combinations involving k^ xt ) Consider the convex combinations 

ai,o,x(x) = Xk^(x) + (1 - X)kf xt (x) = A^ + (1 - A)^, A G [0, 1], 

of the smallest element of JC + and the largest element of JC. Since 

e' /2 ai,o,A(e*) = ^"T* + Acosh(t/2) 
cosn(r/2) 
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is unbounded for any A G (0,1], it cannot be positive definite and hence combining an 
arbitrarily small amount of (the largest element of K~) with the smallest element of K + 
moves out of K + . 



A similar argument can be used to show that any k G K for which the measure m in (2.15 ) 
has the property that m({0}) > cannot be in K + . 

Example 4.3. (Convex combinations of kf xt and A;J; xt ) Replacing /cq x * in the previous ex- 
ample with another k^ xt does sometimes yield convex combination in K + . To be precise, 
let 

a 1)V)X (x) = XK xt (x) + (1 - A) AG [0,1], (4.1) 

1 + x 

be a convex combination of the smallest kf xt (x) = 2/(1 + x) of K + and other extreme points 
k^ xt of /C, v G (0,1). We are interested in the problem to determine A for which ai iU ,x 
belongs to K + . Our result is that for every v G [0, 1), a\^ v ^\ is in K + if and only if 

2 



< A < 



(1 + y^) 2 M/4 + v -l/4) 



(4.2) 



Moreover, a\,v,\ is an extreme point of K + if and only if equality holds in (4.2). 

Since the proofs require some technical results from Section [5j they are postponed to 



6.2 



Note that the right-hand side of ((42j) is < ^ but a 1)1)A = 2/(1 + x) G K + for 

1. 



Section 

A G [0, lj. Thus, this example exhibits some discontinuous behavior at v 

It is straightforward (see the last paragraph of Section 2) to extend these results to show 
that the function 

It 

ai,,,A(x- 1 ) = A/* xt (x) + (1 - A) — - , A G [0, 1), 

1 + x 



is in T if and only if the inequality holds in (4.2) and that it is an extreme point of T if 
and only equality holds. 

Example 4.4. (Extended Heron means) Consider the convex combinations of x~ 1//2 and 
extreme points of /C, i.e., 



Xkl xt + (1 - \)x 



-1/2 



(4.3) 



which are sometimes known as Heron means when v 

/C + for all A G [0,1]. However, for v = the function (4.3) is in fC + on ly for A 



1, in which case (4.3) is obviously in 

since 



x 1//2 /cQ Xt (x) is unbounded. Furthermore, it follows from Proposition 3.7 that for v G (0,1) 
and A / the function in (4.3) is never in /C + because k^ xt K, + and x 1 / 2 k ( ^ t (x) — > as 
x — > oo. 

Next, consider ( |4.3| ) with v = as the convex combination of the largest and the smallest 
elements of K~ for A G (0, ^) . Since 

„2i \ -i 



-t 



1 + e 2 



+ (1 - A)e 



-t 



1 

A 



1 



cosh t + 



l-A 



3.3 



(2), we have AA?f xt + (l-A)x- 1 / 2 K~ 



with (1— A) /A > lis not positive definite by Lemma : 
showing that Kr is not convex. However, the dual set J-~ is convex and / i— >• k = 1/ / 
transforms T~ to JC~ . Thus, although Kr is not convex, harmonic means of functions in Kr 
are in K~ . 
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4.2 Families of classic functions in fC 

Examples given so far suggest us that ZC + is a rather thin subset of /C. Therefore, it is a 
bit surprising that we find a number of one-parameter families in JC + in the examples below. 
Each of these families shows some type of symmetry and monotonicity in the =3! order on 
maximally suitable intervals. In all these cases, the symmetry condition xk{x) = /c(x _1 ) can 
be easily checked and it is rather straightforward to use the Pick mapping property to verify 
that they are in /C. Although the most intriguing family is associated with the WYD skew 
information, it is also rather complex. 

Example 4.5. (Heinz type means) The family of functions 



= x a + X l-a > ^[0,1], (4-4) 



has the dual family 



TH/ , 1 x- a + x~ 1+a 



which were used in [53]. One easily recovers the Heinz type means via (2.10) since 

a G [0,1]. 



y x a y 1 a + x l a y a 



k*(x/y) 2 
In addition to k^ 2 (x) = x~ 1 / 2 = k^ 2 (x), important special cases are 

fc H (x) = kf(x) = y^~ x = kT\x), k$(x) = k?(x) = ±±Z = ft8*(s) 

reflecting the obvious symmetry around x = \. 

Since e t ^ 2 k^(e t ) = 1/ cosh((a — is positive definite, k^ G fC + for any a G [0,1] and 
k& G /C~ for any a G [0,1]. A different proof of the former was in [6} Example 3]. 

If < a < j3 < |, then =4 ( see [13 Section 2]) so that the pair of functions k^ for 
a G [0, |] and k^ for a G [5, l] can be regarded as a single family which increases in the =3! 
order from the smallest to the largest element of JC. 

Moreover, whenever < a < /3 < \, 

feg(e t )_ cosh((i-/3)t) 
fcj(e*) cosh ((I -a)*) 

is infinitely divisible by [8] Theorem 1]. 

Example 4.6. (Binomial means or power means) The functions 

/ 2 \ 1 I° L 
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are easily verified to be in tC as observed in [39, Theorem 3 (i)] and correspond to the binomial 
(or power) means 

1/a 

a €[-1,1]. 



y 



Interesting special cases are 



y 



1 + X 

2x 



feg«(x) 



^ B l/2( X ) 



kft(x) = lim k%(x) = x- 1/2 , 

a— >0 



(1 + Vx) 

4x 
4 



k 



WYD 
1/2 

WYD, 



k l/2( X ) ~ (1 + y^)2 ~ fc l/2 ( X ) 



fcf(x) 



l + X 



kr(x), 



where /c}| VYD is given in Example 4.8 Moreover, k^(x) = k_ a (x) which implies that for this 
family 



a \ — 1 



a € 



-1,1], 



with the obvious abuse of notation. It follows from [27\ Theorem 9] that if— 1 < /3 < a < 1 
then k® =4 kg, so that we have a decreasing family in the =^ order. Since /eg (x) = x -1 / 2 , 
conclude 



we 



• k% G 1C + if and only if a G [0, 1], 

• G /C~ if and only if a G [-1,0]. 

Moreover, k^(e t )/kg(e t ) is infinitely divisible whenever (5 < a \27\ Theorem 9]. 
Example 4.7. (Power difference means) The family of functions 

k P a G (x)^^-- X -^, a €[-1,2], 
a — 1 x" — 1 

give the family of power difference means considered in \19\ 120] . In fact, 

V Ti/r I \_ a — 1 x a — y a 



M a (x,y) 



,,01— 1 ' 



(4.5) 



k^P(x/y) " ' a x a 1 — y L 

whose family is also called the A-L-G interpolation means since it interpolates the arithmetic, 
the logarithmic and the geometric means by allowing us to recover all of these as special cases 

1 + X , on / \ ,. , / \ X — 1 



£?(*) 

I'D 

1/2 



^1/2 ( X ) 



k™{x) 



2x 
2 



fc PD (x) = lim ,^ D (x) = - , 

Q-s>0 x log X 



A- 



PD; 



lim fc PD (x) 

a— >1 



logx 



1' 



l + X 

It is known |X9(, Proposition 4.2] that A; PD G K, for all a G [—1,2]. Moreover, we have 
k^P = kf^ a , which implies that for this family 



Q \ — 1 



1_a 



a G 



•1,2], 



with the obvious abuse of notation. If— l</3<a<2 then fc PD =<! /c? D (see |19t Theorem 



2.1]), so that we have another increasing family. Thus, since & P 5(a 



- 1 / 2 , we can conclude 
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• jfe™ G K + if and only if a G [|,2] , 

• G /C~ if and only if a G [-1, §] . 

Moreover, the monotonicity can be strengthened to the infinite divisibility of £^ D (e*) /fc« D (e*) 
for /3 < a by |27t Theorem 5]. 

Example 4.8. (WYD family) One of the best known families in K, is an outgrowth of the 
Wigner-Yanase-Dyson skew information discussed in Section 2.4 which leads to the functions 

> (1-x) 2 ' P£[ ' J - 1 bj 

This family is symmetric around p = |, and the special cases p = 0, 1 should be understood 
by continuity, i.e. 

uWYD/ n _ ,.WYD/„ n _ i- lWYD/^ 1o § x 



kf u (x) = fc (^) = lim /fc p WYU (x) 



p->l p x - 1 

Other important special cases are 

lWYD/ \ _ 4 lWYD/ \ _ J.WYD / \ _ 4 1 + y/x + X 

*l/2 W-(l + v ^)2' *-l/2W-*3/2 W- 3 ^(1 + ^)2' 

We can summarize the CP situation for this family as follows: 

(a) A;™ G /c + if and only if p G [ 0j 1 j j 

(b) A;™ £ /c - if and only if p e [ _ l, _i] u [f , 2] . 

For p £ the functions k^ YT) increase monotonically with respect to the =^ order. 

Set r = p + q — 1, a = p/r, (3 = q/r so that r > and < a < j3. We note 

fc^ YD (^) = g(i-g) a-x^a-^-p) 

fcWYD^) p(l-p) ' (l-X^l-x 1 ^) 

g(a--l) (1 -x ra )(l -x^" 1 )) _ M a (x r ,l) 
~ a(/3- 1) ' (l-x^Xl-x 1 ^- 1 )) ~ JW^l)' 



where M a (x,y) is the power difference mean defined by (4.5) (for any real parameter a). 
Therefore, when ^ < p < q < 2, fc^ VYD (e')//c^ VYD (e') is infinitely divisible by \27\ Theorem 
5] and in particular /c^ VYD =^ k^ YT) . 

Thus, the functions k^ YT) form a smooth family which are in /C + up to p = 1 when p 
increases from i. Therefore, fe^^ D lies on the boundary of /C + , and we conjecture that it is 
an extreme point of /C + . 

The operator for k = k^ YD is given by (1.1), which implies that /c^ VYD G /C + for 
p = 0,1. In the proof of [6l Theorem 2], explicit (double) integral representations were 
obtained for fi^, when = kp for p G (0, 1) in such a way that the CP of 0,^ immediately 
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follows and hence k^p^ G /C + for p G (0,1). This gives the "if" part of (a). An alternate 
proof of this, as well as details for the remaining claim above are given in Section 6.1 This 
requires results from Section 5.1 which are of independent interest. 

Unlike other families we consider, the functions k^p^fa) = l/fc^ VYD (x _1 ) do not belong 
to the WYD family. Despite the extensive study of WYD metrics, there seems to have been 
little attention given to this dual family 

k™ YD (x)=p(l- P ) 7 r, P€[-1,2]. 

p K 1 FK y> (x - x x ~p)(x - xpy 

This is symmetric around p=\ and special cases are 

By (2.17) and (b) above the functions fc^ YD are in JC+ for p G [|, 2] and in KT for p G [0, 1]. 

Example 4.9. (Stolarsky means) As in the WYD example above, the dual of the Stolarsky 
function gives a different family. In this case, we introduce both 

i . i 

\a(x — 1) J \ a(l — x) J 

It is known [391 Theorem 3 (ill)] (also [6 } Theorem 3]) that fe„ G /C for a G [—2,2] and this 
range for a such that fc-? G /C is optimal. The functions fe^ (x) correspond to the familiar 
family of Stolarsky means as 

ykl t (y/x) = S a (x,y)=[ -. \ ) . (4.7) 



k^(x/y) " ' " ' ' '' V a ( x -y), 

The mean S\(x,y) = e~ 1 (x x /y y ) 1 /( x ~ v * ) for a = 1 is called the identric mean. 
The functions k^" include more familiar special cases than k^ as follows: 



= TT^' kf(x) = limk s a t (x) = ex^, kf /2 {x) = + 4 ^_ 

1/3 



1 + x' lvy Q ^i av; ' ' (1 + v 7 ^) 2 ' 

1 + x N 



fc St (x) = lim kf(x) = ^(x) = A: st 2 (x) = f 

a-+0 X — 1 \ 



2x 2 



which provide an interesting comparison with the other families, as shown in Table 1. 
When 

' (3(x a -y a )\^ 



S a ,/3{x,y) 



a(xP - yP) 



it was proved in [271 Theorem 12] that S a> p(e x , 1)/S a i^'(e x , 1) is infinitely divisible as long 
as a < a' and (3 < (3' . Since in particular S a> i(x, y) = S a (x, y), it follows from (4.7) that this 
implies that when a < (3 the dual family k^ 1 =4 k^ increases and fejg* =4 k^ decreases. We 
can then conclude that 



20 



• ka G K- + and € /C~ if and only if a G [-1, 2], 

• kf G XT and G AC+ if and only if a G [-2, -1]. 

As remarked above, the dual functions form a different family with special cases 

k%(x) = ( * . ) ki\(x) = x- 1 / 2 , kl\x) = lim fc^x) = -^i, 

= (1+ 4 ^ )2 , kf t (x) = hmk s a \x) = e- 1 x^, %!\x) = = k™\ x ). 

The pair kf i i a for — 1 < a < 2 and k^ for —1 < a < 2 can be regarded as a single family 
which increases in the ^ order from k^ xt to kf xt . The functions kf * and /c^ give new members 
of ZC + which do not not appear in any of the other families. Moreover, k^} 2 must lie on the 
boundary of both IC + and /C, which implies that tC + touches the boundary of tC at the interior 
of a face. It seems reasonable to conjecture that k_ 2 is an extreme point of /C + . 



It is interesting to compare the behavior of these examples as the parameter a changes, 
as summarized in Table [T] and Figure [TJ 
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Table 1: Summary of common crossing points 



The Stolarsky family is the only one which goes through all of the indicated points. The 
WYD family is the only one which does not begin and end at the smallest and largest 
elements, and moves outside of both /C + and Kr for some parameter range. 



4.3 Geometric bridges 



In Examples |4.2| and 4.3 of the previous section we considered arithmetic weighted averages 
of 2/(1 + x) (the smallest of /C + ) or x -1 / 2 (the largest of K, + ) with extreme points of /C, and 
noticed that such averages can be in JC + in rather limited cases. In this section we consider 
a different type of averages, often called a geometric bridge, which is defined as weighted 
geometric means [fci(x)] 1_A [A;2(x)] A , < A < 1, of ki,k% G K. We first show that K and 
are all closed under geometric bridge interpolations as far as some infinite divisibility 
condition is satisfied for . The equivalence of (ii) and (iii) in the next theorem implies 
that a similar result holds for T and J-^ . 



21 



Figure 1: Schematic diagram of families in fC parameterized so that they increase in the =^ 
order. The lower ball corresponds to K, + and the upper ball to K,~ . The curves are described 
from right to left, with the first three beginning at the smallest member 2/(1 + x). The 
rightmost curve (red) describes the Heinz family k^ p ; the next (blue) curve the binomial 
family k® p ; the next (green) curve the power difference family ki_ p . The brown curve on the 
left the WYD family kp^ YD in the range p£ [|, 2] and the dotted brown curve on the right 
the dual WYD family. The crossings at 4/(1 + \fx) 2 and logx/(x — 1) can easily be seen. 
The complex Stolarsky family, which is the only one which starts at the smallest 2/(1 + x) 
and goes through both of these crossings while remaining in /C + before reaching x~ 1 / 2 , is not 
shown. 



Proposition 4.10. If k\,k2 G K,, then for every A G [0,1] the function [/ci(x)] 1_ [/^(x)] 
is also in K,. Moreover, if k\,ki G IC + (resp., JC~) and one of the following conditions is 
satisfied, then for every A G [0, 1] the function [/ci(x)] 1 ~ A [fc2(x)] A is also in K, + {resp., Kr) : 

(i) both e*/ 2 fci(e*) and e t / 2 k2(e t ) are infinitely divisible, 

(ii) &2(e*) / ki(e t ) is infinitely divisible, 

(iii) fci(e*)/&2(e*) is infinitely divisible. 



Proof. To prove the first assertion, let k\,k2 G /C; then by Theorem 2.4 they have the Pick 
mapping property, from which it follows that [/ci(x)] 1_A [/c2(x)] A also has this property and 
hence is operator monotone decreasing. Since the symmetry condition in Theorem 2.4 is 
obvious, we conclude that k\~ x k2 G fC. 

To prove the second assertion, let ki,k2 G /C + and < A < 1. When (i) is sat- 
isfied, [e t ^ 2 ki(e t )] 1 ~ x and [e t / 2 k2(e t )] x are positive definite and hence so is the product 
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e t/2 [A;i(e t )] 1 - A [A;2(e')] A . We note 

so that we get the desired positive definiteness from either (ii) or (iii). 



Finally, the assertion for /C is easily verified by taking kj(x) = l/kj(x ) and using 



(2.17). QED 



Recall that all the one-parameter families in tC + given in Examples 4.5-4.9 satisfy the 
property of infinite divisibility (an order stronger than =<!). Therefore, the above proposition 
implies that geometric bridges joining k\, ki in each of these family sits inside KL + . 

Example 4.11. Consider the bridge 

T a/2 _ _— a/2 

ka(x) = [kfix^lx^r = . _ n ■ 



By Proposition 4.10 this is in /C + for a G [0, 1]. In fact, k a G /C in the larger range a G [0, 2]. 
One way to see this is to observe that gp(x) = x~^(l — x) = x~@ — x 1- ' 3 is operator convex 
for /3 G [0, 1]. Then it follows from [331 Theorem 11.13] that 

gp(x) + xgpix- 1 ) xP-x'l 3 

— W^W — = = 2/3 

is a multiple of a function in K, for /3 G [0,1]. 



Since 



/,. , , 2/ ; _ 1 . sinh(at) 
a sinh t 



it is easy to see by Lemma 3.3 (1) that k a G /C + if and only if a G [0, 1] while k a G Kr if 
and only if a G [1, 2]. It is also known [8j Theorem 2] that k a (e t )/kp(e t ) is infinitely divisible 
whenever a < (3. Given the special cases 

ko(x) = JmfcaCar) = *4/ 2 0) = ^ | ^3/4 = fc? /4 C«0, 

fc!(x) = x- 1 / 2 , fc 2 (x) = ^ = fc^x), 

it follows that k a {x) is a family which increases on [0, 2] in the ==! order from \ogx/{x — 1) to 
(l + x)/2x. 



The connection between gp and k 2 p is interesting because, as mentioned in Section 2.4 
g(x) = (x — l) 2 k(x) is always an operator convex function with the properties needed to define 
a symmetric quasi-entropy. Although one can begin with a function g(x) which does not 
satisfy g(x) = xg(x~ 1 ) and generate a function k G JC, it is not at all obvious how to reverse 
the process without obtaining a symmetric g. In this case, we have found an asymmetric 
g, in particular gi/ 2 = x -1 / 2 — x 1//2 , which generates the key function k[x) = x -1 / 2 G K,. 
The associated quasi-entropies do not seem to have been studied previously, but appeared 
recently in 
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The remaining examples are concerned with geometric bridges joining kf xt and other 
extreme points of JC which require a more difficult analysis. 

Example 4.12. For fi, u, A G [0, 1] we define 

with k^ xt given by (2.13). This is in fC by Proposition 4.10 A special case 



was treated in [61 Example 5]. We have 

1-2A 

which is positive definite exactly when < A < \ since 1/ coshi is infinitely divisible (see [HI 
Theorem 1] for instance). Therefore, .91,0, a is in K, + if and only if < A < ^. 

Example 4.13. For the more general case 

91>V>X ( X ) = k?\ X ) ) ={ TT - X ) { {x + v){1 , + vx) ) ■ ( 4 - 8 ) 

which increase pointwise with A G [0, 1] from /cf xt to kf, xt . Its behavior (in the present context) 
when v G (0, 1) seems much more mysterious. Our results here are: 

(i) the pointwise order of g\ )V ,\ m A can be also strengthened to the =^ order, and conse- 
quently the set {A G [0, 1] : gi )V \ G /C + } is a subinterval [0, A c (^)], 

(ii) for each v G (0, 1) the critical value \ c (v) satisfies 

\ < *c(") < \. (4.9) 

The proof requires some lengthy computations of Fourier transforms, which will be presented 
in Section 6.3. Unfortunately, we do not have any information about the form of A c (^). 

Example 4.14. It is worth noting that a family of modified bridges 

*..-*»(«) = + , -a)u + a - Aw )"' 0£A£1 - <" 0) 

joining /cj xt and /cq x * was constructed in [15] for the explicit purpose of finding a one-parameter 
family which increases from kf xt and k^ xt in the pointwise order and all of whose elements 
except for A = 1 are regular (here k G K, is regular if lim^o < +00). Without the 
regularity requirement, the families in Examples 4.5-4.7 and 4.9 have this property in the 
stronger =<! order. 
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Prom the same computation as in Lemma |6.4| (Section 6.3) we observe that the set A 
{ A E [0, 1] : ^ ^- + } includes [0, j] . Koenraad Audenaert did some numerical work 

suggesting that gi t i-x,x is not in /C + for A > 0.3 giving a CP crossing at a point slightly 
smaller than 0.3 which would be consistent with (4.9). However, we do not know strong 
monotonicity in the =^ order for the family ( |4.10 ). To conclude that A is of the form [0, a] we 
would need a stronger result, e.g., that \ c (v) is monotone in v. 



5 Positive definite functions 

In this section, we present results on positive definiteness and infinite divisibility of certain 
functions involving hyperbolic functions, which are needed in our proofs. The study here is 
considered as a continuation of [26} [28l 127] . which are of independent interest. 



5.1 Positive definiteness of sinh ratios 

We investigate positive definiteness of the function 

j,,. _ sinh(af) sinh(6i) 
sinh 2 1 

with a, b > 0. If a, b < 1, then f(t) is a positive definite function as the product of two such 



functions (see Lemma 3.3(1)). It is actually infinitely divisible as is explained in [HJ [26] for 



instance. We will show that the converse also holds true. 

Theorem 5.1. The function f(t) is positive definite if and only if a,b < 1. 

When a + b > 2, we have limt_>±oo f(t) = oo so that f(t) cannot be positive definite. 
When a + b = 2 and a ^ b, the obvious estimate 

/<«)=<*< (2±*)' = !=.]&/(«> 

also shows failure of positive definiteness. 

We will assume a + b < 2, and we must show that f(t) is not positive definite as long as 
a > 1 (and hence < b < 1). For this purpose it suffices to deal with a, b rational. Indeed, if 
fit) were positive definite for such a, b (and the result is known for such rational parameters), 
then with a' , b' rational satisfying 1 < a' < a and < b' < b the product 

sinh(a't) sinh(6't) sinh(a't) sinh(6't) 
^ ' sinh(at) sinh(6t) ~ sinh 2 1 

would be positive definite, a contradiction. 

Hence, we will assume that a, b are rational in the rest. Obviously we can further assume 

a = 2* > 1, b = - > 0, a + b < 2 with n,m,k G N even. (5.1) 



The most delicate part in our proof for Theorem 5.1 is covered in the next lemma, and the 
rest of the subsection will be devoted to its proof. 
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Lemma 5.2. The function f(t) cannot be positive definite for a, b rational described by (5.1 ). 

For a fixed s£lwe set 

F(z) = f(z) e isz (= sinh ( az ) sinh(fez) ^ sz = sinh(^z) sinh (| z) ^ \ 
\ sinh 2 z sinh 2 z j 



with a, b given by (5.1), and compute its integral along the following rectangle V: 

Tx z = t, t:-R^R, 
T2 z = R + is, s : — > nn, 
r,3 z = t + inn, t : R —> —R, 
T4 z = —R + is, s : nn — > 0. 

We observe 

sinh(i + inn) = sinhi, 

sinh (™ (t + inn)) = sinh (~ i + imn) = sinh (™ ij , 
sinh (i + inir)) = sinh (~ i + i/c7r) = sinh (| i) 

(since n, to, A; are even) so that we have F(t + inn) = f(t) e iS (*+ m7r ) an d 

r i-R r-R 

/ F(z)dz= / f{t)e ist dt+ / f{t)e ist e- nns dt 
Jriur 3 j-_r 



(l-e" n ™) / f(t)e ist dt 

J-R 

r R 

2e- n7rs/2 smh{nns/2) / f{t)e lst dt. (5.2) 

J —R 



Since a + 6 < 2, we have /(z) — >• uniformly on the strip {z G C; < Imz < nn} as 
Re z —> ±00 and hence 



/ F(z)dz = 0. (5.3) 



lim 

7r 2 ur 4 

Therefore, the Fourier transform of f(t) can be computed from J r F(z) dz. If the Fourier 
transform fails to be positive, then Lemma |5.2| follows from Bochner's theorem. 

Note that z = 0,inn are zeros of sinh 2 z of order 2. However, these two points are also 
zeros for sinh z), sinh (J^z) so that z = 0, inn are removable singularities for F(z). The 
poles (inside of V) closest to T are 

Z\ = in and z n -\ = i(n — l)n. 

Note that z\,z n _i are zeros for sinh 2 z (appearing in the denominator) of order 2 and that 
they are not zeros for sinh (^ z) and sinh z) (due tol<a = ^<2 and < b = ^ < 1). 
Thus, we conclude that z = z\, z n ^\ are double poles for F(z). 

We begin with computation of the residue Res(-F(z); z\) at z = z\. Thanks to sinhz = 
— sinh(z — in) (or by direct computation) the power series expansion of sinh z around z\ = in 
is given by 

sinhz = -{{z- zi) + (z- zi) 3 /3! + (z - zi) 5 /5! + • • • ) 



26 



= -(z - zi)(l + (z - Zl ) 2 /3l + (z - zi) 4 /5! + •••)■ 
We thus get the following Laurent series expansion: 

1 1 1 



sinh 2 z (z-zi) 2 (i + ( z _ Zl )2/ 3! + ( z _ Zl )4/ 5! + ...) 2 

1 1 

(z — z\) 2 1 + (z — zi) 2 /3 + higher even powers 
(l — (z — ^i) 2 /3 + higher even powers) . 



Since 

we have 
Computations 



dz* 



I Z = Z\ 



(1 + is(z - zt) - s 2 (z - Zl f/2 + •••). 



sinh (™ zij = i sin (™ 7r) (= i sin(ai)) , 



dz 



sinh ( — z) 



^ cosh ( ^ z) 



\z=z\ n 



— COS ( — 7T I 



a cos 



(at)) 



give rise to 
and similarly 



sinh z) = i sin(a7r) + acos(a7r)(z — zi) + 
sinh z) = i sin(&7r) + 6cos(67r)(z — Zi) + • 



From (5.4)-(5.7) the Laurent series expansion of F(z) around z = z\ is given by 
7 t^t (l — (z — zi) 2 /3 + higher even powers) (l + is(z — z\) + • • • ) 

(z — Zl) Z ' y ' 

x (2 sin(a7r) + a cos(a7r)(z — zi) + • • • ) (i sin(&7r) + &cos(67r)(z — z%) + • • • ). 



(5.4) 



(5.5) 



(5.6) 
(5.7) 



(5.8) 



The residue Res(i r (z); z\) is nothing but the coefficient of (z — zi) -1 here, i.e., that of (z — z\) 
in the product of the above four brackets (multiplied by e _7rs ). Since the starting term is 1 
and a (z — zi)-term is absent in the first bracket, what we have to compute is the coefficient 
of (z — z\) in the product of the last three brackets. In this way we arrive at 

Res(F(z); zi) = e~ ns (i sin(a7r) • 6cos(&7r) + acos(a7r) • i sin(&7r) 

+ is ■ i sin(a7r) • i sin(&7r)) 
= ie~ ns (a cos(a7r) sin(&7r) + b cos(6-7r) sin(a7r) — s sin(a7r) sin(i>7r)) . 

We next move to computation of the residue Res(F(z); z n _i) at z = z n _i (= i(n — l)7r). 
Because of sinh z = — sinh(z — (n — X)m) (= — sinh(z + 7T*)) with n even we have 

sinhz = -((z - z n _i) + (z - z n _i) 3 /3! + (z - z n _i) 5 /5! + • • • ) 
= -(z - z n _i)(l + (z - z„_!) 2 /3! + (z - z„_!) 4 /5! + • • • ) 
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(with the identical coefficients as in the power expansion around z = Zi), and hence we have 



sinh 2 z (z - z n -i) 2 



Z n -1, 



2 /3 + higher even powers) 



d Jlsz 



Z=2 n _l 



(is)V 



(isYe-^-V™, (5.5) has to 



again (see (5.4)). Also, since 
be replaced by 

e isz = e -(n-l)« ( X + is{z _ ^ _ s 2 {z _ Zn _ l} 2 /2 + . . . ) 

So far w e ha ve not seen changes of coefficients except the obvious modification that the factor 
e -7rs j n (5 5) was replaced by e ~(™^ 1 ) 7rs . On the other hand, since 

sinh (^ z n -i) = i sin (^ (n — l)7r) = i sin (m7r — ^ 7r) 
= — isin(^7r) (= — isin(ai)), 



sinh (^2) 



^cosh(^z) 



COS m7T 



z=z n -i 

n) = 



III 

n 



^cos(^(n-l)vr) 
cos 7r) (= acos(ai)) 



the power series expansions (5.6) and (5.7) are replaced by 

sinh (™z) = — i sin(o7r) + a cos(a7r)(z — Zn-i) H , 

sinh (| z) = — i sin(&7r) + bcos(b-K)(z — 2 n _i) -\ 

with constant terms of the opposite sign. The four relevant expansions are now at our 
disposal, and the same reasoning as before (see the product (5.8)) gives us the following 
conclusion: 

Res(-F(z); z n -i) = e"*™ -1 ^ 8 (— i sin(a7r) • bcos(bn) — acos(a7r) • isin(&7r) 

+ is (— i sin(a7r))(— i sin(6-7r))) 
"( n_1 ) 7rs (_ acos(a7r) sin(67r) — 6cos(6-7r) sin(a7r) 

— s sin(a7r) sin(&7r)) . 



le 



The sum (multiplied by 2-ki) of the two residues we have computed so far can be rearranged 
in the following way: 

2m(Res(F(z); z x ) + Res(F(z); z n -ij) 
= 2n [e~ ns (—a cos(a7r) sin(&7r) — b cos(&7r) sin(a7r) + s sin(a7r) sin(6-7r)) 



+ e 



-(n—l)ns 



acos(air) sin(&7r) + 6cos(&7r) sin(a7r) + s sin(a7r) sin(&7r))] 



2ne~ nns/2 [ e (n/2 ~ 1)?r ' s (-a cos(a7r) sin(&7r) - b cos(bn) sin(a7r) + s sin(a7r) sin(&7r)) 
acos(a7r) sin(&7r) + 6cos(fr7r) sin(a7r) + s sin(a7r) sin(&7r))] 



+ e 



-(n/2-l)-n\s , 



Aire n7TS / 2 [ — (a cos(a7r) sin(&7r) + b cos(bir) sin(a7r)) sinh ((n/2 — l)-7rs) 

+ s sin(o7r) sin(&7r) cosh ((n/2 — l)irs) 1 . 



Therefore, by recalling (|5.2|) and (|5.3|) we conclude 
1 

2tt" 



1 r 00 
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1 



sinh(n7rs/2) 



sin(a7r) sin(&7r) • s cosh ((n/2 — 1)tvs) 

— (a cos(a7r) sin(&7r) + 6cos(&7r) sin(o7r)) sinh ((n/2 — l)7rs) 
+ lower order terms 



(5.9) 



A few remarks concerning "lower order terms" are in order. Other candidates for poles (inside 
of r) of F(z) are 

Zi = ilir (for£ = 2,3,...,n-2), 

where nature of singularities at these points (i.e., removable singularities or poles of order 
at most 2) is determined according to values of sinh (^2^), sinh (^zA appearing in the 
numerator. Anyway, residues arising from them give us linear combinations of factors of the 
forms 

sinh ((n/2 tts) , cosh ((n/2 - £') vrs) with £' = 2, 3, . . . , n/2 

in the above big bracket ( |5.9[ ) (possibly with the linear factor s for double poles). Indeed, 
the only source for exponential factors is the power series expansions of e lsz around z = z# 



(see (5.5)), which actually gives rise to 

e iszt = e -lns = e -nn S /2 e (n/2-t)ns = 2 , 3, . . . , n - 2) . 



Thus, by recalling the factor e nns / 2 appearing in (5.2), we get the assertion. 
The dominant term (as s — > ±oo) in the numerator of the Fourier transform is 

sin(a7r) sin(&7r) • s cosh ((n/2 — 1)tts) , 

and we observe 

sin(a7r) sin(6-7r) < 



thanks to 1 < a < 2 and < b < 1 (see (5.1)). Consequently, the Fourier transform takes 



negative values for \s\ large (i.e., failure of positive definiteness for /(s)), and Lemma 5.2 has 
been proved. 



5.2 Fourier transform of ((cosh(t/2) + a)(cosht + (3)) 



-i 



Detailed information on positive definiteness for (cosh fc (£/2)(cosh£ + /3) m ) will be needed 
to prove results on geometric bridges in Example |4.13| However, a direct computation for 
its Fourier transform based on residue calculus seems hopeless due to the fact that poles 
of higher orders have to be considered. Instead, in this subsection we compute the Fourier 
transform in the special case k = m = 1 with the additional parameter a as in the theorem 



below (and then in Section 5.4 we will check higher order partial derivatives relative to a and 
(3 to achieve our goal). 

Theorem 5.3. For a £ (—1, 1) and > 1 we have 
1 f°° e ist alt 



4tt y_oo (cosh(£/2) + a) (coshi + f3) 
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1 

sinh(27rs) 



sinh(26>s) 



cr 



(2a 2 -1 + 0) 



@~2^ cos(As) sinh(7rs) — a sin(As) cosh(7rs) 



VW 2 



K^+a 2 ) 



(5.10) 



where 9 = cos 1 a G (0,7r) and A = log (d + ^/ (3 2 — lj , i.e., A > is a solution of cosh A = j3. 
Proof. For a fixed s£iwe set 

p isz 

F(z) = - 

(cosh(z/2) + a) (cosh z + (3) 

and compute its integral along the following rectangle T: 

Ti z = t, t:-R^R, 

T 2 z = R + is, s : -> 47T, 

T 3 z = t + +ni, t:R^r -R, 

T 4 z = -R + is, s : 47r -> 0. 

Due to cosh((i + 47ri)/2) = cosh(t/2), cosh(t + 4-7ri) = coshi and 



lim 

7r 2 ur 4 



/(z) dz = 



we have 



lim / f(z) dz = lim / /( 
r-^oo J T r-^oo Jriur 3 



z) dz 



lim (1 - e" 47rs ) 

R->oo 



R e ist ds 

R (cosh(t/2) + a) (coshi + 0) 



2e- 2ns sinh(2vrs) / 
J— e 



e 4St ds 



(cosh(i/2) + a) (coshi + /3) 



and we will compute J r f(z) dz by residue calculus. 

It is easy to see that we have the following six simple poles inside of T: 

z = 2i{iT-6), zi = 2i(ir + 0), $ = iir + A, £f = 3in + A. 

When a = 1 (i.e., 9 = 0), 2-7ri is a double pole. However, we assumed a 6 (—1,1) to avoid 
this complication. Note that the Fourier transform formula (5.10) itself remains valid for 
a = 1 by the obvious limiting argument with the understanding 



sinh(26»s) 



a=l 



sinh(2#s) 
lim - v ; = 2s 
ayi VI - a 2 



(see (5.17) below). We note 

Res(zj,F(z)) 



2 smh(zj/2) (cosh Zj + 0) 



3 = 0, 1, 
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and observe 



e «s2 _ e -2(n-0)s & isz 1 _ e -2(n+0)s 

sinh(zo/2) = isin(7r — 9) = ism6 = i\J 1 — a 2 , 



sinh(zi/2) = isin(-7r + 9) = —is'm9 = —i\J 1 — a 2 , 

cosh(zo) = cos(2(vr - 9)) = cos(2#) = 2 cos 2 9 - 1 = 2a 2 - 1, 

cosh(zi) = cos(2(vr + 9)) = cos(26») = 2a 2 - 1. 

Thus, we compute 

Res(z ;F(z)) = T 

2 



Res(zi;F(z)) = 



\Vl^^(2a 2 -1 + p) V / T^ 2 (2a 2 - 1 + 0) ' 
e -2(-K+e)s 2ie~ 2( - n+e ^ s 



^Vl^ 2 (2a 2 - l + p) Vl^ 2 (2a 2 - 1 + 0)' 



and consequently we have 

2ie~ 2 ™ (e 26s - e' 2ds ) 



Res(z ; F(z)) + Res(zi; F(z)) 



Vl-a 2 (2a 2 - 1 + 13) 

Aie~ 2ns sinh(2(9s) 
Vl-a 2 (2a 2 + 



We note 



Res(Zf;F(z))= j—^- ■ , . uc± , j = 0,l. 



(cosh(£±/2)+a) sinh^ 

We observe 



cosh(^/2) = cosh((ivr ± A)/2) = ±isinh(A/2) = ±iJ^^ = ±iy 



cosh(£f /2) = cosh((3ivr ± A)/2) = =Rsinh(A/2) = +i 

sinh(^) = sinh(i7r ± A) = =F sinh A = =F\/ /5 2 - 1, 
sinh(^) = sinh(3i7r ± A) = ' P 2 ~ 1, 

and hence 

— 7TS±iAs ■ — 7rs±iAs 

Res(^;F(z))- 
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Res(£f;F(z)) 



-3irs±i\s 



-3ns±i\s 



We compute 



ie 



Res(e + ;F(z))+Res(e -;F(z)) 



»W . //3-1 



+ ia) + e 



ia 



2ie 



y^^cos(As) — asin(As) 



i e 



-37TS 



Re S (e i + ;F(z))+Re S (er;i ? W) 



As/ 



ia) + e 



+ ia 



2ie 



-37TS 



cos(As) + a sin(As) 



x/^T(^ + « 2 ) 

Both the quantities have sin, cos, and we conclude 

Res(£+; F(z)) + Res(^; F(z)) + Res(£+; F(*)) + Res(^; F[z)) 



2i 


^ cos(As) (e" 7 " - e" 37rs ) - a sin 


(As)(e" 7rs + e" 37rs ) 









4ie" 



-2ns 



2 



cos(As) sinh(7rs) — asin(As) coshes) 



(5.12) 



The desired Fourier transform formula (5.10) is obtained as the sum of (5.11) and (5.12) 
(multiplied by 2m). QED 



5.3 Analysis of ((cosh(t/2) + a) (cosh t + (3)) 1 

Here, we recall the Kolmogorov theorem (a version of Levy-Khintchine formula): A function 
f(t) on M is the characteristic function of an infinitely divisible probability measure with 
finite second moment if and only if there exist a finite positive measure v and a 7 £ M such 
that 

, j* / \ r°° fe its - , , s 

log f(t) =vyt+ = dv{s) . 
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Detailed accounts can be found in (37J E] f° r instance. We note that functions f(t) we are 
dealing with here are all smooth and hence the "finite second moment" condition is automatic 
(see (37J §2-3] for instance). 

In the following lemma we state two explicit examples of the Kolmogorov theorem obtained 
in \27\ Lemma 2 (ii) and Lemma 16] for later use and for the convenience of the reader: 

Lemma 5.4. 



(i) For a > and 9 G [0, 7r), 

, 1 + cos ( 

log 



cosh(at) + cos 
(ii) For a > and A > 0, 

log 



1 + cosh A 



Ats 



cosh(at) + cosh A / 
When a = (i.e., 9 = vr/2), ( |5.10[ ) reduces to 



\ cosh(#s/a) 

ist) — ——. j—r- as. 

ssmh(7rs/a) 



s cos(As/a) 

ist) — as. 

ssmh{irs/a) 



e ist dt 



cosh(t/2) (cosht + /3) 



2vr 



L - y / ^tCos(As) 
(/3 — 1) cosh(-7rs) 



(>0), 



(5.13) 



which corresponds to the special case a = in the next result ([HI Theorem 4.13] and see 
also [271 Section 7]). 



Corollary 5.5. We set 



G(t) 



(cosh(t/2) + a) (cosht + 0) 
with a, P > —1. 

(i) When (3 > 1, G(t) is positive definite if and only if a G ( — 1,0]. 

(ii) When —1 < f3 < 1, G(t) is infinitely divisible for each a G (— l,oo). 



Proof. Assume (3 > 1. Due to (5.13) G(t) is positive definite for a = and remains so for 
a G (—1,0] as well thanks to positive definiteness of 

cosh(i/2) _ —a 



cosh(t/2) + a 



1 + 



cosh(t/2) + a 



(see Lemma |3^3|(2)). When a G (0,1), we have 9 = cos" 1 o_G (0, vr/2) in ( |5.10[ ). Thus, 
the dominant terms in the big bracket in the right side of ( 5.10[ ) are cos(As) sinh(7rs) and 
sin(As) cosh(7rs) so that the quantity in the big bracket takes both positive and negative 
values for \s\ large. To prove (i), it remains to show failure of positive definiteness for a > 1. 
However this follows from positive definiteness of 



cosh(t/2) + a 
cosh(t/2) + \ 



1 + 



a 



cosh(i/2) + \ 
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for instance (and the already known failure of positive deflniteness for a = i). 

Next, assume —1 < (3 < 1. The statement (ii) is obvious for a G (—1, 1], G(t) being the 
product of two infinitely divisible functions under these circumstances. On the other hand, 
when a > 1, we have 



log((l + a)(l + 0)G(t)) 



Ast 



. i cos(29s) cosh(As) , , 

1 - ist) ■ w 9 \ + ■ A \ I ds 
1 ssmn(27rsj ssmh(7rsj 



with 9 = log ya + \/a 2 — lj and A = cos 1 /3 (by Lemma 
written as 



cos(26>s) + 2 cosh(7rs) cosh(As^ 



s sinh(27rs) 



5.4). The density here can be 

(5.14) 



which is certainly positive. QED 



5.4 Analysis of (cosh fc (t/2)(cosht + /3) m ) 1 



In this section we obtain a result which will be used in Theorem 16.61 of Section [6.31 to obtain 
a bound for the interval in which g± u \ is in /C + . 



We assume j3 > 1 and a G (—1,0] as in Corollary 5.5 (i). Under these circumstances the 

cosh(26's) + 2 cos(As) cosh(-7rs) 



density (5.14) is switched to 



log((l + a)(l + /3)G(t)) 



1 



ist) 



ds 



s sinh(2-7rs) 

with 9 = cos -1 a G [vr/2, ir) and A = log (/3 + y 7 /? 2 — 1 j . Thus, the positive definite function 



G(t) (Corollary 5.5 (i)) is infinitely divisible if and only if 

cosh(2#s) + 2 cos(As) cosh(7rs) > 0, 



s G 



This is quite a delicate condition, but for the extreme value a = (i.e., 9 = tt/2) the condition 
simply means 

(1 + 2 cos(As)) cosIi(tts) > 0, 
and it is never fulfilled for any (3 > 1 (which is exactly \27\ Theorem 15]). 

We take higher order partial derivatives d^~ 1 d^ c ~ 1 from the Fourier transform formula 



(5.10) (with the variable s fixed). It is obvious that from the left side we get a scalar multiple 
of 

e ist ds 



(cosh(i /2) + a) k (cosh t + 0) ' 
and for the special value a = the above integral reduces to 

e ist ds 

, cosh fc (£/2)(coshi + /3) m ' 



(5.15) 



Therefore, behavior on the Fourier transform (5.15) can be seen by computing d 1 ! 1 1 of 



the right side of (5.10) at first and then by substituting a = 
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The right side R(a, p) of the formula ( 5.10| ) consists of the three terms: 



„n „ / n \ sinh(2#s) „ , „. . sinh(7rs) „ . „. . cosh(-7rs) 

R(a, 0) = F (a, p) . - . v o ; - F c (a, 0) cos(As) . - / 7 . + F s (a, /3) sin(As) 



sinh(2-7rs) 



sinh(27rs) 



sinh(2-7rs) 



_ , „. sinh(20s) „ , „. cos(As) _ . sin(As) 



sinh(27rs) 



2 cosh(7rs) 



2 sinh(7rs) 



(5.16) 



with 



F Q {a,0) 



VT^{2a 2 -1 + 0)' 



'P-1 



n 



+ a 2 ) ' 



We note 



da 



V^l(^ + a 2 ) 



dA 



sin ( 



(5.17) 



and consequently 

9 a sinh(26>s) 

sin(As) = 



2scosh(2(9s) 



5 Q cosh(2#s) 



2ssinh(20s) 



s cos(As) 
7^ 



1 



cos(As) 



s sin(As) 

7W~ 



i 



We begin with the first term in (5.16). Since sinh(2#s) and cosh(2#s) behave like constants 



against dp, d^ t ~ 1 d^~ 1 of the first term is a polynomial of s of degree at most k — 1 with 
coefficients sinh(2#,s), cosh(2(9s), 1/ sinh(27rs) and so on. Therefore, the substitution a = 
(i.e., 9 = iK I '2) gives rise to a polynomial of s of degree at most k — 1 with coefficients 
containing 



sinh(26»s) 



sinh(27rs) 



cosh(26>s) 



=71-/2 



2cosh(7rs)' sinh(2vTs) 



1 



e=w/2 



2 sinh(-7rs) 



The same procedure for the second and third terms in (5.16) obviously gives rise to a 



polynomial of s of degree at most m — 1. It is important to make sure that the order is 
exactly m — 1, and we will closely check the coefficient of s m_1 . For this purpose we begin 



with the third term in (|5.16|) and we note 
F s (a,p) 



1 



1 



2y^ 



1 



a + i 



+ 



a — i 



(3-1 



d k - x F s (a,P) 



+ 



a + 1 



B-l 



a — i 



B-l 
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d k a - l [F s (a,P) 



sin(As) 
2 sinh(7rs^ 



sin(As) 
2 sinh(7rs) 



So far no s-terms show up because A just depends on /3. We then take derivatives relative to 
p. It is plain to see that the highest s m_1 -term arises from 



1 / sin(As) 



2 sinh(-7rs) 



d k a - l F s {a,P) 



dj- l sin(\s) 
2 sinh(7rs) 



From (5.17) we also easily observe 

m—l 



d^' 1 sin(As) = I 



I (/3 2 -l)( 
By substituting a = 0, we observe 



^77^ \l TTaT sin(As) + lower s-terms (for m odd), 

(RZ — ]\(m— 1)/2 

±tt; ttttt cos (As) + lower s-terms (for m even). 

2 n(m-l)/2 



(5.18) 



a=0 



V 



,9-1 



+ 



/3-1 



^0 



/ 



as long as A; is even. From the discussion so far, for k even the highest s m 1 -term arising 
from 



,»<-). ;)/.- 1 / ... sin(As) 



\ 2 smh(7rs) 

is a non-zero scalar (of course depending upon 0) multiple of 



a=0 



„m— 1 



sin(As) 



sinh(-/rs) 
depending upon the parity of m 



(for m odd) or 



„m— 1 



cos(As) 



sinh(7rs) 



(for m even) 



We next move to the second term in (5.16). We note 



2 a l M 2 



0-1 



( 



d k a - l F c (a,(3) 



(_l)fc-i(A;-i)! 



V 



a + i 



/3-1 



a — 1 



0-1 



J 



The presence of the minus sign this time in the big bracket enables us to conclude 



d k a - l F c (a,f3) 



a=0 



^0 
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for k odd. Since the formula akin to (5.18) is available to cos(As), for k odd the highest 
s m -term arising from 



\ 2 COSh(7TS) 



a=0 



is a non-zero scalar multiple of 



„m— 1 



cos(As) 



,m-l 



COsh(7Ts) 



(for m odd) or 



sin(As) 



cosh(-7rs) 



(for m even) 



this time. 



Summing up the discussions so far, we conclude: For s large the lead ing t erms of 
d^ l ~ 1 d^ c ~ 1 R(a, (3) (which arise form the first term and the last two terms in (5.16) respec- 
tively) are non-zero scalar multiples of 



„fc-i 



and 



s™" 1 sin(As + 6) 



e"° e" 

regardless of the parity of k. 

We are now ready to prove the following result: 

Theorem 5.6. We assume j3 > 1 and set 

H(t) ^ ' 



cosh fc (£/2)(cosht + /3y 
for positive integers k,m. 

(i) H(t) is positive definite if and only if k > m. 

(ii) H{t) is infinitely divisible if and only if k > 2m. 

Proof, (i) Firstly we assume k > m. Since 

1 1 / 



cosh fc (t/2)(cosht + /3) m cosh fc - m (t/2) I cosh(f/2)(cosht + /?) 



positive definiteness of H (t) follows from Corollary 5.5 (i) (or rather the paragraph before the 



corollary). On the other hand, when k < m, for s large the dominant term in (5.15) (i.e., the 
Fourier transform of H(t)) is 



,m— 1 



sin(As + 5) 



as was mentioned in the paragraph right before the theorem. Thus, the Fourier transform 
admits both positive and negative values and hence H{t) cannot be positive definite. 



5.3 



By Lemma 



(ii) We set A = log [ft + ^/ f3 2 - lj > as in Theorem 

/•CO 

log((l + /3) m H(t)) = / (e ist - 1 - ist) F(s) ds 



5.4 



we have 
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with 

k mcos(As) k + 2m cos(As) 

F(s) = ~t~ = . 

2ssinh(7rs) ssinh(7rs) 2s sinh(7rs) 

Thus, H(t) is infinitely divisible if and only if k + 2mcos(As) > 0, i.e., > 2m. QED 

Note that the optimal case in (ii) (i.e., case k = 2m) corresponds to infinite divisibility 
of ((cosht + l)(cosht + (see [271 §7.1]) because of cosh 2 (i/2) = (cosht + l)/2. The 

function l/cosh(i/2) is positive definite (and indeed infinitely divisible) while l/(coshi + (3) 
(with f3 > 1) is not (see Lemma |3.3| (2)). Thus, intuition might suggest that as far as the 
function H{t) is concerned one has higher (resp., lower) chance for positive definiteness and/or 
infinite divisibility as k (resp., m) increases. The theorem completely clarifies where proper 
balance is taken. 



6 Proofs of results in Section [4] 



6.1 Results on WYD family 



For the functions fc}^ YD , p£ [—1,2], defined by (4.6), we will prove the next results stated in 



Example 4.8 



Theorem 6.1. 

(a) The function fc^ YD belongs to /C + if and only if ' p G [0, 1], 

(b) The function fc^ VYD belongs to K.~ if and only if p G [ — 1,— |] U [|,2] . 

Proof. We may and do assume p£ 2] in view of the symmetry /c^ VYD = /cJY_ YD , and set 

, U \ - *j,wyd/ 2ix _ 1 sinh(pt)sinh((l - p)t) 
f p {t)-ek p (e )- p(1 _ p y s[Qh 2 t 

g p (t) = et/k^(e-*)=p(l-p). ^ ' 



sinh(pt) sinh((l — p)t) 



for convenience. Theorem 3.4 says that we have to determine when f p (t) and g p (t) are positive 
definite. 

Positive definiteness for f p (t) with p G [3,1] is well-known (where fi(t) is understood as 
tj sinht). When p G (1,2], we have 

1 sinh(pi) sinh((p — l)t) 



which fails to be positive definite thanks to Theorem 5.1 (due to the presence of sinh(pt) with 
p > 1 in the numerator). Thus, we have shown that f p (t) is positive definite if and only if 
p G [g, l] , that is, (a) is proved. 

To prove (b), we next check positive definiteness for g p (t). When p G l] , we have 

lim g p (t) = +00 

t— >±co 
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(where g\ (i) is understood as sinh i/i) so that g p (t) cannot be positive definite. We now move 
to the case p£ (1,2] so that we use the expression 

sinh 2 t 

9p{t) = P(P ~ 1) 



sinh (pi) sinh((p — l)i) 

When 2 > p + (p — 1) (i.e., p < |), g p (t) once again diverges as i — )• ±oo and fails to be 
positive definite. Thus, it remains to show positive definiteness for p > |. For the extreme 
value p = I we compute 

3 sinh 2 i 3 sinh 2 i 

53/2 (*) 



4 sinh(3i/2) sinh(i/2) 4 sinh 2 (i/2) (4cosh 2 (i/2) - l) 
3cosh 2 (i/2) 3(coshi + l)/2 3 



T 1 + 



4cosh 2 (i/2) - 1 2(coshi + l)-l 4 I coshi + 1 



Since 1/ (coshi + 1) is positive definite (see Lemma 3.3 (2)), so is 53/2 (£)■ Finally, the obvious 
identity 

sinh 2 i sinh 2 i sinh(3i/2) sin(i/2) 



sinh(pi) sinh((p — l)i) sinh(3i/2) sin(i/2) sinh(pi) sin((p — l)i) 

3 . 

2'- 



takes care of the remaining case (i.e., p G (|,2l). QED 



6.2 Proofs for Example 4.3 



We now prove the claims in Example 4.3 



Theorem 6.2. We assume v G (0, 1) and A G [0, 1]. Then, the function the function ax jV \ 

(6.1) 



defined in (4.1) belongs to K. + if and only if 

2^ 



A < 



(l + ^f (1/1/4 + 2 ' 



Moreover, a\ jU: \ is an extreme point in K, + if and only if equality holds in (6.1). 
Proof. With /3 = (1 + v 2 )/2v (> 1) we compute 



e i/2 a!.,,A(e*) 



wfl 1 n cosh(i/2) 1 

A(/? + 1} coshi + /3 + (1 " A) cosh(i/2) 

A(/3 + 1) (coshi + 1) /2 + (1 - A) (coshi + /3) 

cosh(i/2) (coshi + /3) 
(A(/3 - l)/2 + 1) (coshi + 0) - A(/3 2 - l)/2 

cosh(i/2) (coshi + /3) 
A(/3-l)/2 + l A(/3 2 -l)/2 



cosh(i/2) 



cosh(i/2) (coshi + j3) 
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Let us recall (5.13) (where the symbol A for cosh 1 /3 there is changed to a to avoid the 
obvious confusion) and 

r°° e ist dt _ 2vr 
J_ 00 cosh(i/2) cosh(7rs)' 
The Fourier transform is thus given by 



1 f°° t/2 ,tMst* A(/3-l) + 2 A(/3 + l)(l-^ m cos(a S ) 
- / e i / 2 ai^ jA (e t ) e ist dt = -± — —— ^ -j — r 

7T J^oc COSh(7TSJ COSh(7TSj 



N(s) 



COsh(7Ts) 

with the numerator 



N(s) = 2(1 - A) + A v / 2(/3 + l)cos(as). 
This computation says that e t / 2 a\ iVj xifi t ) is positive definite (i.e., ai,v,\ £ ^ + ) if and only if 

Av/2(/3 + l) <2(1-A). 

It is obviously satisfied for A = (which corresponds to the obvious positive definiteness of 
e t//2 fc| xt (e*)) while for A > the requirement is the same as 



Finally, we compute 



l + /3 +1= + + 



2 



2 V Au ' 

which proves the first part. 

Now, we set A(z^) = 2-^/(1 + \Jv) 2 and prove the second part. It is obvious that a\^^\ is 
not an extreme point of /C + if A < A(z^). To prove the converse we assume that ai ^xtu) ( x ) = 
Xki(x) + (1 — X)k2(x) with some A 6 (0, 1), ki, ki G fC + and hence 

ki(x)= k e u xt (x)dmi(v), i = l,2, 
J[o,i] 

with representing probability measures mi on [0, 1]. From the uniqueness of a representing 
measure we have 

Ami + (1 - A)m 2 = (1 - \{v))b\ + \{y)b v . 

In particular, we have 

Am a ({l}) + (1 - A)m 2 ({l}) = 1 - \{y\ 

Ami ({»/}) + (1 - A)m 2 ({^}) = X(y), ^ J 
and they sum up to 

A(mi({l}) + mi(M) + (1 - A)(m 2 ({l}) + m 2 ({i/})) = 1. (6.3) 
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On the other hand, we have mj({l}) + mi({f}) < 1 (i = 1, 2) because m^s are probability 
measures. Therefore, (6.3) guarantees mj({l}) + mj({^}) = 1 (i = 1,2). Hence, both of 
mi,m2 are supported on the two-point set {u, 1}, that is, of the form 



mi 



(1 - di)5i + 



1,2) 



with dj = 777j({^}) G [0, 1]. Since fcj G /C + (« = 1, 2), the first part of the theorem implies ctj < 
A(V) and hence the second equation of (6.2) forces a\ = 02 = A(V), i.e., k\ = ki = ai iV ,\(v)- 
QED 

Since 2y/u/ (1 + v^) 2 < | (with z/ 7^ 1), we have «i,i/,a ^ + as long as A > i (and 

ve[Q,l)). Choose the extreme value A = \{v) = (^/(J3 + l)/2 + l) * with /? = (1 + t; 2 )/2ia 
Then, it is straightforward to compute 



/3 + 1 cosht+ v / 2(/3 + 1) - 1 



2 cosh(t/2) (cosht + /3) 
It is not clear if this positive definite function is infinitely divisible. 



6.3 Results on geometric bridges 

We consider the geometric bridges gi^ : \ (x), < A < 1, between fe| xt and kf, xt with v G [0, 1) 



given by (4.8) and described in Example 4.13. Since the case v = was settled in Example 



4.12:, we assume v G (0, 1). Our main result is 



Theorem 6.3. For each fixed v G (0, 1), there is a critical A c (dependent onv) such that the 
function gi :Ut \ is in K, + for A G [0, A c ] . Moreover, 



^ < A c (z/) < - for each v G (0, 1). 



(6.4) 



This will follow from a series of lemmas and theorems below, which are of independent 



interest. First, observe that one can rewrite (4.8) as 

2 



9i,u,x(x) = x l ' 2 — 



1/2 + 3.-1/2 



1-2A 



(l + vf 



2v [ + 



2v 



and define for f3 



2v 



> 1, 



1 



cosh 1_2A (i/2) Vcosht + /3 



Recall that Theorem 3.4 implies that g\ v \ G /C + if and only if f u ,x(t) is positive definite. 



Lemma 6.4. The function fu,\(t) is infinitely divisible if and only if < A < \. 
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Proof. With a = log (^(3 + a//3 2 — 1^ > by Lemma 



5.4 



we have 



log f V) \(t) = A log 



1 + 
cosh i + /3 



+ (1 - 2A) log 



cosh(t/2) 



(< 



V.S'f 



isi) F(s) ds 



with 



, , cos(as) , ts 1 
F(s) = A . I \ + (1 - 2A) 



2A(cos(as) - 1) + 1 



ssinh(7rs) 2ssinh(7rs) 2ssinh(7rs) 

The minimum of cos(as) — 1 is —2 so that the above density F(s) is non-negative exactly 
when -4A + 1 > 0. QED 

For fixed v the functions gi, u ,\ increase monotonically with A in the =^ order. 
Lemma 6.5. If X' < X, then gi iU> \' =4 9l,v,\ 

Proof. It suffices to show positive definiteness of 

fv,\> jt) _ gi^X'je 1 ) 
fu,\(t) 9i,u,x( et ) 

for A' < A. However, this ratio is equal to 

A— A' 



(6.5) 



1 



cosh 2 ( A - A ')(t/2) V 1 + /? 



cosh t + j3 



1 



cosh t + j3 
cosh 2 (t/2) ' 1 + 13 



2 cosht + /3 
1 + j3 ' cosh t + 1 



A-A' 



A-A' 



On the other hand, by Lemma 5.4 (ii) we have 
2 cosht + /3 N 



log 



1 + cosh t + 1 



\ 1 — cos(as) 

1 — ist) — — r- as 

' s smh(7rs) 



with a = log \J ft 2 — 1 ) > 0. The density (1 — cos(as))/s sinh(7rs) here is positive so 



that (cosht + /3)/(coshi + 1) is infinitely divisible and hence the ratio (6.5) (with A' < A) is 
positive definite. QED 

The monotonicity shown above implies that for each fixed v, the set 

{AG [0,1] : 0i,„, A e/C + } 



is a subinterval [0,A C ], for which we prove that the critical value A c = X c {y) satisfies (6.4) 
above. The lower bound j follows immediately from Lemma 6.4 The upper bound follows 
immediately from Theorem 6.6 below. 

Theorem 6.6. When A > |, the function gi >L/ ,x does not belong to /C + for any v G (0, 1). 
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Proof. By Theorem 3.4, this is equivalent to showing that f u ,\(t) is not positive definite when 
A > 3, for which we will prove by contradiction. We choose a rational ^ with 



1 m 

o < - < A. 
3 n 



(6.6) 



Then by Lemma 6.5, positive definiteness of f u ,\{t) implies that so is Since f Ut m(t) 

is equal to 



cosh 1- ~^(i/2) (cosh t + /3) " 
up to a positive constant, its nth power 

1 

cosh fc (t/2)(cosht + /3) m 

with k = n — 2m would be also positive definite. However, this contradicts Theorem |5.6| (i) 
because of (6.6), i.e., k < m. QED 

How the critical value A c (i/) depends on v G (0, 1) seems to be an interesting problem. We 
note that f Vj i/z(t) is equal to 



cosh 1 / 3 (t/2) fcosht+ to 



1/3 



up to a positive scalar. Although the function (cosh (t/2) (cosh t + (1 + i/ 2 )/2z^)) 1 is known 
not to be infinitely divisible ( |271 Theorem 15] and see also Theorem 5.6 (ii)), its cubic root 
might be positive definite for some values of u). However, the authors are unable to handle 
this delicate phenomenon. 



Acknowledgments 

This work was begun when FH, DP and MBR were participating in the Fall, 2010 program 
at the Mittag-Leffler Institute in Sweden. They are grateful for the opportunity and support 
provided there. The work of MBR was partially supported by NSF grant CCF 1018401. 
FH and HK acknowledge support by Grant-in-Aid for Scientific Research (C)21540208 and 
(C)23540215, respectively. We thank Ben Criger for providing the figure. 



A Proofs from Section [2] 



A.l Proof of Proposition 2.2 



The following is from the proof of \22\ Corollary 2.2] . Since k is an operator convex function on 
(0, oo), the function g(x) = k(l+x) is operator convex on (—1, 1). By Kraus' theorem [29] (see 
[21 Lemma III. 1] or [71 Theorem V. 3. 10]) the divided difference function h(x) = (g(x)—g(0))/x 
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is operator monotone on (—1, 1). Then by Lowner's integral representation (see [21 Theorem 
II. 1] or [71 Corollary V.4.5]) there exists a (unique) finite measure [i on [—1, 1] such that 



x 

h(x) = a+ I — dfi(X) 

[-1,11 1 ~ Xx 



with a = h(0) = fc'(l). Thus 

f x 2 

(1 + x)k(l + x) = (1 + x)(ax + b) + (l + x) / - — d//(A) 

7[-i,i] 1 - Ax 

= (1 + x)(ax + b) + m({-1}) + x 2 / -Hf - dAt(A) 

with b = k(l). Since (x + 1)(1 — Ax) < 1 whenever x G (—1, 1) and A G (—1, 1], the Lebesgue 
dominated convergence theorem implies that 



lim / ^±f d M (A) = 
z\-l J i-iA] 1 - Ax 



x-\-i 7r_i i] 1 — Ax 

so that 



lim x k{x) = lim (1 + x) k(l + x) = //({ — 1}). 

x\0 a:\-l 



A.2 Proof of Theorem EOl 

Since the divided difference function /i(x) = (A;(x) — k(l))/(x — 1) is operator monotone on 
(0, oo) as in Section A.l it is known (see, e.g., [13, Theorem 1.10]) that there exist a (unique) 
7 > and a (unique) positive measure fi on [0, oo) with Jj 0oo \(l + A) -1 dfi(X) < +oo such 
that 

x — 1 



h(x) = h(l) + 7(x - 1) + / rdfi(X), xG(0,oo). 

'[0,oo) z + A 



k(x) = /t(l) + fc / (l)(x-l)+ 7 (x-l) 2 + / ^ t^^(A), xG(0,oo). (A.l) 

J[0,oo) x + A 

the symmetry property xk(x) = /c(x _1 ) we notice that 

/cf x) f x 2 
lim = lim x _2 /c(x _1 ) = lim x 2 k(x) = lim / dfi(X) = 

:r->-oo x x-+oc x 7 x\0 x\0 Jin X + X 



by the Lebesgue convergence theorem. On the other hand, since (x — l) 2 /x(x + A) / 1 as 
x /* oo, we have 

lim 1 / (X - \) d/x(A) = / dji(A) 



^oox7[ 0jOo) x + A y [0jOo) 
by the monotone convergence theorem, and hence 

k{x 



lim 

x— >-oo x 



k'(l) + 7 • (+oo) + f dfi(X). 

J[0,oa) 
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Therefore, 7 = 0, fj, is a finite measure, and k'(l) + Jj » dfi(X) = 0. From 

(^~l) 2 _ x 1 (*-!)(! + A) 

x + A x + A 

it follows that 

k(x) = k(l)- [ {X ~ ^ + A) tfyt(A) 

J [0,oo) X + A 

= ^)+ / f^-l)(l + A)d/i(A), 
J[o,oo)V2; + a / 

which is operator monotone decreasing since so is (1 + A)/(x + A). Moreover, since (x — 
1)(1 + A)/(x + A) / 1 + A as x /■ 00, we have < fe(l) - / [0oo) (l + A) d/i(A) so that 
i|o 00) + ^) d[J,(\) < fc(l) < +00. Now, defining a finite positive measure f on [0, 00] by 

di/(A) = (1 + A) dp(A) on [0,oo), u({oo}) = k(l) - f (l + A)d/x(A), 

J[0,oo) 

we write 

k(x)= [ i^dz/(A), i£(0,oo), (A.2) 
■/[O.oo] x + A 

where (1 + X)/(x + A) = 1 for A = 00. Letting du(X) = dv{X~ l ) on [0, 00] we also write 
k(x) = f du{X) = [ dm 

J[Q,oo] x + A 7[ ,oo] 1 + Ax 

so that 

*(*-') = / 

•/ [0,oo] X + A 

This is the familiar integral expression of the operator monotone function fc(x _1 ) with a 



unique representing measure v. Hence the measure v satisfying (A.2) is unique (this fact 
itself is also well-known). Since 

it follows that v = v. Define 

dm{X) = 2dv{X) on [0,1), m({l}) = u({l}), 

to obtain 

t( *H,>(^ + ^H + ^ ({i}) 



l l + x (! + A ) 2 j 

y M («+»)(!+».) ■^ _ * B(A) - 
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Finally, note that the uniqueness of m is immediate from that of v in (A.l). 



The integral expression (A.l) was also given in 
|13| Theorem 5.1]. 



which was considerably extended 

It 



m \16\ ineorem b.lj. There is another route to prove the two theorems in Section 2.1 
was proved in [U Theorem 3.1] that a function k : (0,oo) — > (0,oo) is operator monotone 
decreasing if and only if it is operator convex and non-increasing in the numerical sense. 
It is easy to see that if an operator convex function k satisfies the symmetry condition 
xk(x) = &(x _1 ), then it is non-increasing numerically. Hence we have the implication (a) 
(b) of Theorem |2.4[ All other parts of Theorem 2.4 are plain or well-known. Then we 



can prove Theorem |2.3| by applying the familiar integral expression to a symmetric operator 
monotone function /c(x _1 ) as above (indeed, this part of the proof is the same as the proof 
of [311 Theorem 4.4]). 



B Contraction bounds 



As stated in (2.7), monotone Riemannian metrics contract under the action of quantum 
channels, i.e., CPT (CP and trace-preserving) maps. It is well-known |21| |4"4"1 [33], I55j that 
the quasi-entropies H g (A,B) = H g (A, B, I) given by (2.8) with K = I contract under CPT 
maps, i.e., 

H g {${A)MB))<H g (A,B), A,B<E F d , 
whenever g is operator convex on (0, oo). In the rest of this subsection let g(x) = (l — x) 2 k(x) 



with k G K, as in Section [2~4| In applications, the maximal contraction rate plays an important 
role, which motivated in [33] the following definitions of contraction coefficients : 

,r Ent W -sup sup H dmm 

d P ,jev d ,p^ Hg(p,-y) 

and 

Vk(&) = SUP SUp SUp rk(Y y-N • 

A contraction coefficient was also defined |51] for the trace norm ||X||i = Tr \X\ = Tr (X*X) 1 / 2 
distance which also contracts under CPT maps, i.e., 

_Dob/*^ „„„ II^(P-7)I|1 



r ? JJoD ($) = sup sup 

d p, 7 e© d ,p^ 7 IIP ~7lli 

where the superscript reflects the fact that this is the quantum analogue of the classical 
Dobrushin coefficient. 

For any CPT map it was shown in |33l Theorem IV. 2] that 

v f [em ($) < ?^ clEnt ($) < 1 for any k € JC, 
and in \5i\ Theorems 13, 14] that 



C i% nt ($) < v Uob (^) < V^ icm ( <i> ) ( B - 1 ) 
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when k(x) is given by (4.4). The upper bound in (B.l), given in |51| Theo rem 3] for the 
particular case kQ Xt (x) = (1 + x)/2x and in [53] for k = k^ in Example 4.5 holds for any 



k G JC. Our work here was motivated by the lower bound in (B.l) based on the following 



observations from |33|, [54"] . Applying the max-min principle to the eigenvalue problem 



( $ n* w 



*)(x) = \n k Jx) 



(B.2) 



(for which X = I always yields the largest eigenvalue Ai = 1) implies that 



r ? ^($) = sup sup A£(*,p), 

d p£V d 



where $ is the adjoint of $ (with respect to the Hilbert-Schmidt inner product) and X^i^, p) 
denotes the second largest eigenvalue of ( |B.2 ). This is equivalent to the eigenvalue problem 



t*,(#(jo) = (tfr 1 o $ o 0* = AX 



for the trace-preserving map $ = (0*) o $ o ^ restricted on H^. When $ is 
positivity-preserving, 



\%($,p) = sup 



lixih 



< sup 



\\x\u 



Dob, 



A sufficient condition for ^ to be positivity-preserving is that both (Qp)' 
CP0 which we have seen holds if and only if k(x 



-1/2 



There may be particular maps 



<]? for which is positivity-preserving even when Qh and/or its inverse are not. Whether 



or not the bound rj (<&) 



< r] b (&) holds for other k G /C even though T*$ is not 
positivity-preserving is an open question. 



C Some pedestrian arguments 

In this section we present, for the benefit of non-experts, some very pedestrian ways to see 
certain well-known results used in this paper. 

C.l Functional calculus for L D and R D 

It is basic that when D has the spectral decomposition D = J2j w j\£.j)(£j\ (where 
is the physicists's notation for the spectral projection onto the eigenspace of the eigenvector 
tp(D) = ^^j) I £?')(£?' I f° r anv function ip(x) on (0, oo). However, writing the spectral 
decomposition of Ld is a bit awkward. This can be avoided by observing 

ip(L D )(A) = J2<P(™Mi)(t;M and ^(R D )(A) = ^(<ML|£,}<0|. 

k k 



5 Unfortunately, in [53] it was claimed that is positivity-preserving for k(x) — logx/(x — 1). Although 
ftp given by (1.1 1 is clearly positivity-preserving, the inverse (fi*) 1 given by (1.2 1 is not. 
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Then the product 

{^{r d Ml d ){x) = 5>(u>i)^K-) 16X01 (&,x$j) 

id 

Since Ld and Rd commute, it follows that for an arbitrary function cj)(x,y) 

4>(L D ,R D )(X) =J2<t>(wi,Wj) i&Atj) 
hi 

which is exactly the Hadamard product of A o X when dij = 4>(u)i, Wj) and X is represented 
in the basis 

C.2 Integral representation and inversion of BKM operator 



Although it is well-known that and its inverse for k(x) = logx/(x — 1) are given by (1.1) 



and (1.2), most proofs rely on an expansion in eigenvalues. Therefore, we give an explicit 



argument here. First we use the elementary formula 



r '°° ( 1 1 
log x = I I du 







1 + U X + u 



to write 



{\ogL D R- D l ){X) = (log Ld — log Rjj) (X) = {L logD -R logD )(X) 

oo / 1 1 \ 

X X dt 

D + tl D + tl J 

oo 1 1 

DTT^ DX - XD) DTi I dt 



from which it follows that 



nUx) = log ^- logi V) = r -J— x -J— dt. 

m ' L D -R D K ' J D + tl D + tl 



Now, observe that 



ii v nl-t dt 



y Z? < [(log£»)X-XlogD]L> 1 " < (it = J j^XD 1 

= DX-XD = (L D - R D ){X) 



so that L>< n k D (X)D 1 - t dt = X, which implies (PHI, i.e., 



= C rfYD^dt. 
Jo 
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